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ABSTRACT 

We have compiled a catalog of optically-selected quasars with simultaneous observa- 
tions in UV/optical and X-ray bands by the Swift Gamma Ray Burst Explorer. Objects 
in this catalog are identified by matching the Swift pointings with the Sloan Digital Sky 
Survey Data Release 5 quasar catalog. The final catalog contains 843 objects, among 
which 637 have both UVOT and XRT observations and 354 of which are detected by both 
instruments. The overall X-ray detection rate is ~ 60% which rises to ~ 85% among 
sources with at least 10 ks of XRT exposure time. We construct the time-averaged 
spectral energy distribution for each of the 354 quasars using UVOT photometric mea- 
surements and XRT spectra. From model fits to these SEDs, we find that the big blue 
bump contributes about '-~^ 0.3 dex to the quasar luminosity. We re-visit the aox-L 
relation by selecting a clean sample with only type 1 radio-quiet quasars; the dispersion 
of this relation is reduced by at least 15% compared to studies that use non-simultaneous 
UV/optical and X-ray data. We only found a weak correlation between iboi/iEdd and 
auv- We do not find significant correlations between a^ and aox, ckox and auv, and c^x 
and log L(0. 3-10 keV). The correlations between auv and a^, aox and ax, aox and auv, 
iboi/iEdd and ax, and Lbd/iEdd and aox are stronger amongst low-redshift quasars, 
indicating that these correlations are likely driven by the changes of SED shape with 
accretion state. 

Subject tieadmgs: quasars: general — catalogs — ultraviolet emission — X-rays 
1. Introduction 
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Variability is a ubiquitous phenomenon 
of quasa rs (e.g., Matthews & Sandage 196 



Smith fc Hoffieit 1963; GaskeU fc Khmek 200 
Vanden Berk et all 120041: IWilhite et all 12005 . 



20061 : iMeusinger et al.l 120111 ) and has been 
observed in radio, infrared, UV/optical, X - 



ray and even 7-ray bands (e.g.. | Kovalev et al 

lnnnollr>: P_ Al :„,„llnnnol Io^i.^-l. „-i- _i 



12002^; 'Rieger & Aharonian"2008':'Sakata et id_ 
i2010i: iPark et al.i 2010; Grupe et al.i i201^ 



While variability provides considerable infor- 
mation on the size scales of the quasar cen- 
tral engine and can be utilized to estimate 
the mass of the central super m assive black 
hole (e.g., IPeterson et al.l 120041 ). it is a sig- 
nificant source of scatter in multi- wavelength 
correlations. Clearly, when fluxes of variable 
sources measured in different wavebands at 
different times are combined, the combined 
spectral shape may not be representative of 
the spectral shape at a specific time. It is 
therefore essential to take simultaneous ob- 
servations when determining multi- waveband 
properties in order to understand the true 
spectral shape. 

The dispersion in multi-waveband corre- 
lations produced by variability can be esti- 
mated. For example, there is evidence that 
at least 60% of the dispersion of the Bald- 
win Effect and at least 75% of the EW(C iv)- 
Qfox relation can be attributed to variability 
(|Wu et al.l 120091 ). However, a large multi- 



waveband program of simultaneous observa- 
tions to verify that the dispersions of various 
relations can be reduced has not yet been per- 
formed. In this work, we will use simultane- 
ously observed UV/optical and X-ray data to 
study quasar SEDs. 

The Swift Gamma-ray Burst (GRB) Ex- 
plorer ( Gehrels et al.l2004[ ) includes co- aligned 
X-ray and UV/optical detectors. The X- 
ray telescope (XRT; Burrows et al. 2005) 
is an imaging spectrometer that covers the 
0.3-10 keV band. The Ultraviolet Optical 
Telescope (UVOT; Roming et al. 2005) pro- 
vides photometry in six bands from ^ 1928 A 
to - 5468 A. By default, XRT and UVOT 
are operated simultaneously. Because Swift 
is primarily a GRB mission, there are few 
quasars that were specifically targeted for 
observations. However, due to the relatively 
large fields of view (FOVs) of both the UVOT 
and the XRT, a large number of serendip- 
itous sources, including quasars, are ob- 
served. The Sloan Digital Sky Survey (SDSS; 
York et al. 2000) Data Release 5 (DR5; 



Adel man-McCarthy et al. 2007) quasar cat- 
alog (jSchneider et al.l l2007t ) contains 77,429 
optically-selected quasars in 5,740 deg^; the 
UVOT FOV is 17' X 17', so there is a density 
of approximately one SDSS quasar per UVOT 
field. We have matched the Swift pointings 
from launch to June 2008 {^ 3.5 years) with 
the SDSS DR5 quasar catalog, to examine the 
properties of quasars observed simultaneously 
in UV and X-ray light. 

Catalogs usin g similar st ra tegies were 



constructed 



iGrupe et al. 



tnici) 



Tueller et al.l (|2008D and 



Our work is unique in 
terms of the sel ection criteria and sa mple size. 
The sample in iTueller et al.l ( 20081 ) contains 
153 hard X-ray (14-195 keV) selected local 
AGNs with a mean rcdshift of ^ 0.03. Be- 
cause hard X-ray photons more easily pen- 
etrate gas and dust than UV/optical and 
soft X-ray photons, this study provides a 
more homogeneous sample of quasars than 
those based on UV/optical a nd soft X-ray 



bands . Using this sample, iWinter et al 
(|2009l ) estimated the fraction of "Hidden" 



AGNs (|Ueda et al.ll2007l ) in the l ocal universe 
to be ~ 24%. The sample by ICrupe et al 



(|2010t ) contains 92 soft X-ray selected AGNs 
with redshifts ranging from 0.002 to 0.349. 
Due to their selection criteria, objects in this 
sample are X-ray bright Type 1 AGNs with 
few Seyfert 1.5 objects (e.g., Mkn 841; see 
Wilkes et al. 1999). Their selection criteria 
are not biased for or against RL AGNs (about 
10% of AGNs are RL). Using this sample, 
Grupe et al.l (|2010l ) constructed composite 



AGN SEDs with simultaneous observations 
from Swift. By fitting these SEDs with two 
different models, they attempted to constrain 
AGN bolometric corrections (BCs). They 
found a significant correlation between UV 
and X-ray spectral indices, namely auv and 
ajj for AGNs with ax > —1.6. Because their 
sample contains a large fraction of narrow line 



^The spectral indices are defined as /,/ oc 1/+" in this 
work. 



Seyfert 1 (NLSl) galaxies, they were also able 
to examine differences between broad line 
Seyfert 1 (BLSl) galaxies and NLSl galaxies 
in terms of the auv-iboi/^Edd relation. 

As a result of the detection limits in their 
selected wavebands, both of the above stud- 
ies are limited to nearby AGNs and rela- 
tively small sample sizes. We find 1034 SDSS 
quasars within 20' of Swift pointings, which 
is almost an order of magnitude larger than 
either of these samples. 

One of the prominent features in the AGN 
SED is the big blue bump (BBB) in the ex- 
treme UV (EUV) energy band, which is be- 
lieved to be primarily produced by thermal 
emiss ion from an accretion disk (e.g., [Shields 
119781) . Because of strong Galactic and (pos- 
sibly) intrinsic extinction, it is almost impos- 
sible to observe this wavelength region. This 
feature is likely, however, to be an important 
contributor to quasar bolometric luminosity. 
Based on the observed data in UV/optical 
and X-ray bands, we can place some con- 
strain ts on the flux contr ibution of this fea- 
ture. iGrupe et al.l (|2010l ) were able to con- 
strain the BBB using their sample of ^ 100 
AGNs at low redshift z < 0.4. In contrast, 
over 50% of the quasars in our sample are at 
z > 1, thus, we will be able to better constrain 
the BBB feature. 

It should be emphasized that we are not 



tral hardne ss between the UV and X-ray 



repeat ing t he global SED work o f lElvis et al 



(|1994[) and [Richards et al.l (|2006l ). We are at 
tempting to use simultaneously acquired data 
to constrain the BBB. This component, be- 
cause of the lack of observational data, is usu- 
ally represented as a power-law connecting 
a UV flux point, such as 25 00 A, to an X- 



ray p oint, such as 2 keV fe.g.. lRichards et al 
20061 1. The slope determi ned by these two 
point s is defined as Oox (|Tananbaum et al. 



"ox = 0.3838 log [i2 kev/-^25oo aJ ^^^ 
which is used to characterize the spec- 



bands (e.g.. lAvni fcTananbaumI 11982 



Anderson &: M argon'l987 t I Wilkes e t al."l994: 



1986 



Vigiiah et al.. 2003: Strat cva et al.ll2 005: Stc ffen et al. 
2006l : IJust et al.ll2007t ). The expres sion above 



follows the notatioin convention in Just et al 



( 2007 ). so L ! is the monochromatic lu- 

^ "' 2500 A 

minosity at 2500 A in ergs s~^ Hz~^. A 
long-standing problem with the measurement 
of ttox is that it is difficult to obtain simulta- 
neous measurements of an object in the X-ray 
and UV/optical bands. Source variability in- 
troduces scatter into measurements of aox- 
Our simultaneous observations remove this 
noise. 

This paper is organized as follows. In Sec- 
tion [21 we describe the observations and data 
processing; in Section [31 we present our data 
processing results, including the UV/optical 
light curves, and the composite SED for each 
quasar; in Section[4l we present our final cata- 
log; in Section[5l we select a sample of quasars 
from this catalog to study the dispersion of 
relation and other correla- 



2500 A, 



the tto 

tions. In Section[6l we summarize this catalog 
and our conclusions. Throughout this work, 
we adopt the following cosmology: flu = 0.3, 
f^A = 0.7, Ho = 70 km s~^ Mpc"!. 

2. Observations and Data Processing 

Our quasar sample was compiled in the fol- 
lowing steps. 

1. Candidate objects for our catalog were 
selected as any SDSS DR5 quasar that 
lie within 20' of the center of the 
Swift FOV in any pointing from launch 
through June 2008. 

2. XRT data were processed to obtain X- 
ray count rates, spectra and spectral pa- 
rameters. 

3. UVOT data were processed to obtain 
UV and optical photometry. 



4. UVOT photometry were supplemented 
with measurements at other wave- 
lengths from published catalogs. 

5. Quasar SEDs were constructed. 

6. Additional parameters were calculated 
based on the SEDs of each quasar. 

The raw sample is constructed by match- 
ing 3.5 years Swift pointings and the SDSS 
DR5 quasar catalog and contains 1034 ob- 
jects. Fig. [1] shows the distribution of our 
quasar sample in the luminosity (represented 
by absolute magnitude Mi) - redshift dia- 
gram. The distribution of our sample in this 
diagram is consistent with the SDSS DR5 
quasar catalog. 

2.1. XRT Data 

We processed the XRT data using the 
task xrtpip eline (HEADAS version 6.10) from 
FTOOLS (|Blackburnlll995[) . For each obser- 
vation identified by an observational ID (here- 
after OBSID), this routine stacks all available 
XRT snapshots (individual exposures with 
the same segment number) observed in pho- 
ton counting mode and generates a composite 
sky image for each OBSID with an associated 
exposure map. 

Source and background extraction regions 
are then defined for each quasar in the candi- 
date list. The quasar coordinates are adopted 
from SDSS DR5 quasar catalog, which are ac- 
curate to « 1". The source and background 
regions are circles with default radii of 30" 
and 120", respectively. These circles are vi- 
sually inspected in the deepest sky image. 
This ensures that the background does not in- 
clude serendipitous X-ray sources. The posi- 
tion and size of the background region circle is 
adjusted so that it lies entirely within the field 
of view in all sky images covering this object. 
We exclude sky images in which the source 
region is not fully within the FOV or not cov- 
ered by all individual exposures. We also dis- 
card sky images in which the source is severely 



contaminated by nearby X-ray sources. De- 
pending on the brightness of the X-ray source, 
the source region radius is adjusted to be large 
enough to ensure that the photon density at 
the border is at the same level as the back- 
ground. We then enlarge the background re- 
gion so that its radius is at least four times as 
large as the source region radius. 

XRT data are prepared using standard 
FTOOLS packages. Event files for each OB- 
SID are extracted and cleaned using XSe- 
lect (version V2. 4a). xrtmkarf (version 0.5.6) 
is used to create an auxiliary response file 
(ARE) for each OBSID. The ARE file in- 
cludes corrections for filter transmission, vi- 
gnetting, effective area and point spread func- 
tion (PSF). It also accounts for hot pixels 
and hot columns, which are masked out from 
data and decrease the effective exposure time. 
The response matrix file (RMF) we use is 
swxpc0tol2s0.20010101v011.rmj^ 

The total number of background-subtracted 
X-ray counts, A'xph, is determined for each 
quasar by combining the counts measured 
in each sky image. This value is used 
to assgin a quality fiag to each object: g 
(good with A'xph > 100), a (acceptable with 
10 < Axph < 100), w (weak with A^xph < 10), 
and (out of FOV) . In the candidate sample 
of 1034 quasars, 103 objects are flagged as 
good {g), 296 are flagged as acceptable (a) 



^We noticed that this RMF file was intended for use 
with data taken from launch through 28 February 
2007 only and that other RMFs are appropriate for 
data taken from March through August 2007 and for 
data taken since 30 August 2007. To test the sensi- 
tivity of our simple spectral models (described below) 
to the choice of response matrix, we fitted models to 
several high signal to noise sources using each RMF 
and found the impact upon the model parameters to 
be small: the spectral index changes by less than 0.02 
and the intrinsic column density changes by no more 
than a few tens of percent. We therefore use a sin- 
gle RMF for simplicity and we choose the one for data 
taken prior to 28 February 2007 as the majority of our 
data is from this time interval and most of our sample 
members were observered at least in part before this 
date. 



and 406 are flagged as weak (w). The rest 
(229 objects) are flagged as o because they 
lack useful sky images. 

Next, we produce summed event lists for 
the source and background regions, including 
all observations. From these we extract time- 
averaged spectra. We also generate the ARF 
file for the composite source spectra using ad- 
darf (version 1.2.6). The X-ray energy spec- 
tra are binned using FTOOLS task grppha 
(version 3.0.1). 

The X-ray energy spectra are fit using 
XSPEC (Arnaud 1996; version 12.5.1n). The 
binning strategy and the statistical method 
we use are described in Table [T] Basically, we 
use x^ statistics and have at least ten spec- 
tral bins if iVxph> 100. If A^xph< 100 we ap- 
ply minimal binning to eliminate unoccupied 
spectral bins (this avoids having bins with 
negative backgroui id-subtract ed fluxes) and 
use Cash statistics (|Cashlll979t ). If A'^xph< 10, 
we only estimate flux or flux upper limits. For 
each object with A^xph^ 10, we fit each X-ray 
spectrum with an absorbed power-law model. 
This model includes Galactic column density 
along the line of sight to the quasar TVh.g and 
possibly additional absorption intrinsic to it 
TVhj. The values of TVh.g are fixed based upon 
the Leiden/Arge ntine/Bonn (LAB) survey 
o f Galactic H l teartma nn k Burton! Il997t 



Kalberla et al.M2005 ). We define four varia- 
tions of this model, depending upon whether 
the photon spectral index (F) or A^n.i are al- 
lowed to vary as a fixed parameter. We deflne 
four variations of this model, depending upon 
whether this spectral index a^ or iVjj.i are 
allowed to vary (models A-D; Tabled). We 
visually inspect the fitting quality and apply 
the following rules to select the model that 
produces the best constraint and quality: 

1. If TVxph < 30, we use Model A. In 
this case, it is impossible to constrain 
the intrinsic absorption and the con- 
straint on dx is poor (Icr uncertainty 
of ttx is larger than ~ 0.5). There- 



fore, we fit only the fiux levels and fix 
a^— —1, which is the average value 
of X-ray spectral index for RQ type 1 
quasars (e .g ., iN andra & Pounds 1994^ 
120031 



Page et al 
20091) . 



I2OO3L I2OO4I : lYoung et all 



2. If 30 < iVxph < 100, we allow one addi- 
tional free parameter, either iVH,i or ay_, 
but not both: 

• If the intrinsic column density iVn,! 
is inconsistent with zero and well 
constrained, we choose Model D. 

• If iVjj.i is consistent with zero, and 
Qfox is reasonably well constrained 
(l^ttxl < 0.5), we choose Model B. 

• If neither of these criteria are met, 
we choose Model A. 

3. In the case that A^xph > 100, if a^ 
is well constrained and the la lower 
bound of TVh.i is not consistent with 
zero, we use Model C. Otherwise, we 
follow Rule 2. 

The visual inspection process is performed by 
three people and a consensus is reached to en- 
sure objectivity. The number of quasars se- 
lected for each fitting model is shown in Ta- 
ble [5] Examples of the fitting are shown in 
Fig.H 

Based on the spectral fitting results, we 
calculate a number of parameters, including 
the photon count rate, the observed and un- 
absorbed X-ray flux between 0.3 and 10 keV, 
and the monochromatic luminosity at 2 keV 
in the emitted frame. These values are in- 
cluded in the final catalog (see Section |4]). 

2.2. Weak X-ray Sources 

We process weak X-ray sources (flagged as 
w in the final catalog) with TVxph < 10 sepa- 
rately because some are not significantly de- 
tected by XRT and we can only estimate their 
flux upper limits. 



We apply the Bayesian method of ([Kraft et ah 



19911) to determine whether an X-ray source 
is detected or not. We define a source to be 
a non-detection if 3a confidence level lower 
limit is consistent with zero source counts. 98 
out of 406 weak X-ray sources are detected 
by XRT. 

For the detected sources, we fit Model A 
to determine F(0.3-10 keV) and /(2 keV), 
which is the obsersved flux between 0.3 and 
10 keV and the observed flux density at 2 keV. 
For undetected sources, we consider two 
cases. If the background-subtracted count 
rate is positive, we can use the XSPEC to 
apply Model A. This has the advantage of 
including the calibrations present in the ARF 
file, but the model will be fitted to the mea- 
sured counts and not the upper limit and 
therefore must be rescaled. The flux upper 
limit is determined using 



F, 



uplim 



F, 



N, 



Xspcc ■ 



uplim 



A^xph 



(2) 



in which A^upiim is the 3cr source photon count 
upper limit. 

If an undetected X-ray source has too few 
counts (typically A^xph< 2), XSPEC is un- 
able to apply models to the data. In these 
cases we manually apply two corrections to 
convert the count limits into limiting count 
rates. The first correction is for vignetting. 
Vignetting is a reduction of the effective area 
of the telescope at off-axis positions. We 
adopt a vignetting function of y(0) = l—CB"^, 
where 6 is the off-axis angle in arcminute. The 
coefficient C varies with energy (Cusumano & 
the XRT Calibraqtion Team, 2006; Kennea, 
private communication); here we adopt the 
value appropriate for 2 kcV as it is representa- 
tive and the uncertainty in our measurements 
is dominated by small number statistics. This 
correction is applied by reducing the effec- 
tive exposure time by the vignetting factor, 
thereby increasing the count rate limit. The 
second correction is to account for the finite 
size of the aperture used to measure source 



counts. We generate an ARF file appropri- 
ate for the size and position of the source ex- 
traction region. The output of this process 
includes a report of the fraction of the source 
fluence enclosed within the region; we divide 
the count limit by this fraction to apply the 
PSF correction. Finally, the vignetting- and 
PSF-corrected count rate limit is converted to 
limit ing observed fiux by using the tool pimms 
(jMu kai 1993) together with the assumption of 
Model A. From this model we determine the 
absorption-corrected flux, flux density and 
monochromatic luminosity at 2 keV. 

2.3. UVOT Data 

Instead of using a pipeline to obtain a co- 
added event file as was done with the XRT 
data, we process each individual UVOT sky 
image for each object. The composite pho- 
tometry is obtained by summing over photon 
counts from each individual image and divid- 
ing the sum by the total exposure time. 

First, we must identify suspect and de- 
fective sky images. Because we are using 
serendipitous observations, it is inevitable 
that some sources are too close to the edge of 
the UVOT FOVH. The goal of this process is 
to exclude the low quality sky images so that 
the default source and background regions 
can be directly applied on good ones. 

The source region is defined using a cir- 
cle with a radius of 3", which is recom- 
mended in the UVO T photometric calibration 
(jPoole et al.ll2008l ). Since a 5" radius aper- 
ture, which contains 85.8 ± 3.8% of the PSF, 
was used for calibrating the UVOT, an aper- 
ture correction is applied to the data when 
running uvotsource (Poole et al. 200 8). The 
inner and outer radii of the background re- 
gions are r = 27"5 and R = 35"0 respectively; 
this is the standard background region used 



•^Objects are not necessarily always within the FOV be- 
cause the matching radius used to select the raw sam- 
ple is larger than UVOT FOV and/or because UVOT 
was working under other modes. 



to construct the firs t GRB afterglow catalog 
( Roming et al.l2009() . The large difference be- 
tween the source and background radii ensure 
that the background area is at least 50 times 
as large as the source region to provide an 
accurate background subtraction. 

We then flag each sky image based on the 
position of the object in the image frame and 
the aspect keyword value. We only use images 
whose aspect values equal DIRECT. The other 
images (flagged as —1) do not have correct as- 
pect corrections and may lead to inaccurate 
photometric results. The detailed image flag- 
ging strategy is tabulated in Table [31 We only 
accept UVOT images flagged with 0, 1 and 2. 

Due to relatively small photon counts or 
positioning uncertainties, a fraction of object 
images are not symmetric or well aligned with 
the default source region circle with a radius 
of 3". It is necessary to co- align the quasar 
image with the center of the source region 
circle because the aperture correction process 
assumes a symmetric photon loss outside the 
pre-defined aperture. To co-align the loca- 
tions, we extract a region of 11 x 11 pixels 
(about 5" X 5") centered on the SDSS quasar 
coordinates. If the total photon count in- 
side this region is less than 20, we regard 
this image as "faint" and assign a fiag of "1" . 
Even if a single image is regarded as "faint" , 
it may contribute to the total object counts 
when stacked with other images (e.g., Lehmer 
et al. 2007). If an image is not fiagged as 
"faint", we attempt to fit the source with a 
two-dimensional Gaussian profile the centroid 
of which is set to be free. The background 
is determined by averaging the photons in- 
side the region centered at the SDSS coordi- 
nate and bounded by two squares with sizes 
of 21 X 21 pixels and 11x11 pixels. If the pho- 
ton count inside the 11x11 pixel region is less 
than three times that of the background level, 
this image is also treated as "faint" (fiagged 
as "1"). Otherwise, we fit the 11 x 11 pixel 
region with a 2-D Gaussian profile. The new 
Gaussian centroid is compared with the origi- 



nal SDSS quasar coordinate. We denote their 
angular separation as S. 

There are three categories of images. 

• If (5 < 0".618, the two coordinates are 
consistent and the default source region 
circle is used. This image is flagged as 
"0". 

• If 0".618 <5 < 3", this image is flagged 
as "2" and the source region circle is re- 
centered at the Gaussian centroid. 

• If (5 > 3", this image is fiagged as "3". 
This large offset could be caused by 
the unrecorded aspect problem, which 
requires manual correction or non- 
Gaussian photon distribution. Less 
than 1% of sky images are flagged as 
"3" , which are dropped. 

In the next step, we deflne the default 
source and background regions for each ob- 
ject, and select the sky image with the longest 
passband central wavelength and the longest 
exposure, because, stellar contamination from 
host galaxies is larger in optical than in UV. 
If the V band is not available or its expo- 
sure is too short (< 100 seconds), we exam- 
ine the B band. If the B is unavailable, we 
check the U band, and then successively the 
UVWl, UVM2 and UVW2 bands. For each 
object, we attempt to exclude all stars, galax- 
ies and bad pixels in the background annulus. 
If the source region is significantly contami- 
nated or most of the background region must 
be masked, we remove this object from the 
UVOT analysis. 

In most cases, the procedures above are 
sufficient to define a source and background 
region for all sky images. However, there re- 
main a number of observations that require 
visual inspection. It is sometimes necessary 
to define a different background region to ex- 
clude the sky images in order to obtain ac- 
ceptable photometry. For example, a number 
of images contain stellar ghost rings. Due to 



differences of aspect and exposure time, the 
influence of these factors may vary for each 
sky image even for a given object. In ad- 
dition, although the ASPCORR keyword may 
be set to DIRECT, the image may stiU suffer 
from uncorrected aspect problems so that a 
single point source may have multiple images 
or even smeared images. Consequently, it is 
necessary to visually inspect each sky image 
to perform the second order region customiza- 
tion and image selection. After the visual 
inspection process (> 50,000 images), 3183 
(~ 1%) UVOT sky images were excluded from 
our UVOT image set. 

We use uvotsource (HEADAS version 6.10) 
to calculate photometry on individual sky im- 
ages using the curves of growth from the Swift 
CALDB for the aperture correction model. 
The composite photometry uses the photon 
counts of all available images. For a given 
waveband, the mean magnitude is calculated 
as 

(TO) = Zpt-2.51og(i?Lss) (3) 

in which Z^t is the zero-point magnitude (the 
magnitude when the count rate is 1 pho- 
ton s~^). The values are taken from Table 6 

is the source 
aperture. 



by multiplying the LSS photon count rate by 
the conversion factors in Table 10 (conver- 



of Poole et al 



mm- -Rlss 

with coincidence-loss, 



count rate 

and large-scale sensitivity corrections applied. 
The associated uncertainty in this magnitude 
measurement is 



(4) 

(^LSS) - y^^T (5) 




and 



^(Alss) 



I]('^iAfl,Lss_J 



Er. 



(6) 



where Ti is the exposure time in image i. We 
also calculate the flux density at the effec- 
tive wavelength for each fllter for each object 



sion f rom power-law spectra) of iPoole et al. 
(l2008h . 



The above equations can be applied to an 
individual sky image as well as a group of sky 
images. This process produces three measure- 
ments for each fllter: photometry for each in- 
dividual sky image, for all sky images in an 
OBSID, and for all OBSIDs of a fiher. 

A total of 675 objects in our sample are 
detected by UVOT. 

3. Quasar SEDs 

We classify objects in the database into 
four types based on data availability (see Ta- 
bleS]). In the description below, "useful data" 
includes both detections and non-detections. 

• Type A contains 637 objects that have 
useful data from both UVOT and XRT. 
Among them, 345 are detected by 
UVOT and XRT. 

• Type B contains the 38 objects that 
have useful data from UVOT but not 
XRT. These objects cannot be included 
in correlation analysis. 

• Type C contains 168 objects that have 
only useful data from XRT but not 
UVOT. We will only create the XRT en- 
ergy spectra for these objects. 

• Type D contains 191 objects that are 
bereft of any useful data from UVOT 
or XRT. They are not included in the 
flnal catalog. 

3.1. Supplementary Data 

As the object redshift increases, the UVOT 
wavebands are progressively shifted into the 
EUV band shortward of Lya; as a conse- 
quence, UVOT data for some quasars cannot 
be used to measure the rest-frame UV power- 
law. The SDSS and 2MASS data, observed 



at longer wavelengths, can be used to extend 
the available measurements to UV and op- 
tical bands in the quasars' rest-frames. We 
supplement our UV photometry with the five 
bands (ugriz) from the SDSS DR5 quasar cat- 
alog and thre e bands ( J, H, JtTs) of the 2M ASS 
(|Cohen et al. 2003; S krutskie et al.ll2006l ) sur- 
vey, when available. In the best cases, we have 
14 photometric data points in UV/Optical 
wavebands: six from Swift/IJVOT, five from 
SDSS, and three from 2MASS. 

3.2. Initial SED Plots 

The monochromatic luminosity corre- 
sponding to each UVOT filter can be eas- 
ily computed after shifting the flux density 
to the quasar's rest-frame. The frequency 
band widths are calculated by converting rest- 
frame FWHMs of corresponding wavebands. 

The SDSS quasar catalog provides mag- 
nitudes that we must convert to flux densi- 
ties. Instead of the asinh magnitude used in 
the general SPSS ph otometric measurements 
(JLupton et al.l 119991 ) . we convert the SDSS 
flux densities from corresponding ba nd mag- 
nitud e using the Pogson definition (jPogson 
18571 ) which has a much simpler analyti- 
cal expression nii = — 2.51og/i//o, where 
/o is the zero-magnitude flux. The Pogson 
magnitude system deviates from the asinh 
system for faint objects. According to the 
asinh softening parameters (h coefficients) 
table in the SDSS DR5 photometric cali- 
bration documenlcl, the difference between 
Pogson and asinh magnitudes is less than 
1% for objects brighter than g — 22.60 and 
u = 22.12. The majority of our objects 
are much brighter than the limit. The dif- 
ference between the Pogson and asinh sys- 
tems is therefore negligible. There is one 
object (SDSSJ122740.85-I-440604.7) whose B 
band magnitude is fainter than 22.60 and 
two objects (SDSSJ020316.37-074832.1, and 
SDSSJ133613.62-K025703.8) whose U band 



magnitudes are fainter than 22.12. They are 
all at high redshift (z > 2.5) with extremely 
low X-ray photon counts (A'^xph < 10), so 
they do not play significant roles statisti- 
cally in the correlation analysis result. For 
wavebands in SDSS, we ad opt the central 
wavele ngth and FWHM from iFukugita et al. 
()l996l ). For wavebands in 2MASS, we adopt 



the isophot al central wavelengt hs and band- 
widths from lCohen et al.l (|2003() . 

For X-ray sources with less than 100 pho- 
ton counts, we use a different binning strategy 
in XSPEC from the one presented in Table [1] 
to avoid large error bars. We use the com- 
mand setplot rebin to re-bin the spectrum 
until each bin has a detection at least as large 
as 1(7 and no more than 100 bins may be so 
combined. Error types are set to quad which 
sums in quadrature the errors on the origi- 
nal bins. The rebin command only affects 
the plot appearance but not the fitting re- 
sults (see Fig. [5] for examples) , but produces 
a clear representation of the average flux lev- 
els and associated uncertainties. Because of 
the relatively large uncertainties in the X-ray 
energy spectrum, we calculate the lower and 
upper error bars of each flux point, instead of 
applying standard error propagation. 

The initial quasar SEDs are generated 
based on all available UV/optical photomet- 
ric data points and X-ray energy spectra. 
Before shifting into the rest-frame, we apply 
Galactic reddening corrections to all avail- 
able wave band flux densities using the stan- 
dard E{B — V) d ependent extinction curve 
(FitzpatricyilMS')- Valu es of E(B - V) ar e 
calculated following Schlegel et al.l ( 19981 ). 
Some examples of quasar SEDs are presented 
in Fig. H 

These initial SEDs cannot be used for SED 
fitting for two main reasons. 



• Broad emission lines can contribute sig- 
nificantly to broad band filter measure- 

ments, which may lead to incorrect SED 

Ihttp : //www . sdss ■ org/dr5/aigorithms/f luxcai . htmi#asiiih.tabT^ shapes. For example, the band covering 



Lya can alter the UV slope by up to 0.2. 

• The dates of the observation from 
SDSS, 2MASS, and Swift for a given 
object differ, often by several years in 
the rest-frame. The data may need to 
be shifted to mitigate the effect of vari- 
ability. 

3.3. Emission Line Correction 

Emission line corrections can be performed 
by subtracting the broad line contribution 
base d on their av erage equivalent width (EW) 
(e.g.. lElvial2012l in preparation). Because we 
have SDSS UV/optical spectra for all quasars, 
we use a more sophisticated method in which 
we convolve the response function R{X) of 
each filter with the observed spectrum /tot (A) 
and the power-law only spectrum /pi (A) to 
calculate the emission line correction factor. 



EC = log 






pi 



(7) 



in which 



i^tot = / /tot(A)i?(A)dA 

Fp, = J fp,{X)R{X)dX 

Because the quasar UV power-law usu- 
ally extends from Lya to ^ 5600 A (e.g., 
Vanden Berk et al.l 120011 ). we only need to 
perform emission corrections to filters cov- 
ered within this wavelength range. However, 
low-redshift quasar spectra do not cover Lya 
and high-redshift quasar spectra do not cover 
5600 A. In these cases, the emission line cor- 
rections are perfo rmed based on the compos- 
ite spectrum in IVanden Berk et all (12001). 
Although the shape of the composite spec- 
trum may not be exactly the same as the real 
spectra, the mean corrections are sufficient 
for power-law slope estimation. 

To estimate the errors introduced by per- 
forming emission line corrections using the 



composite spectra, we examine the emission 
line correction trends as a function of red- 
shift for different filters and find that the 
emission line correction is typically less than 
0.1 dex (see Fig. S]). The real spectrum cor- 
rections are generally distributed around the 
composite spectrum, meaning they are gen- 
erally in agreement. To view their differ- 
ences more clearly, we plot the distributions 
of differences between these to corrections, 
AEC = ECrcai - ECcomposito, in Fig. M (only 
the SDSS g band and UVOT B band are dis- 
played, but other filters are similar). These 
distributions indicate that on average, the 
composite spectrum emission line correction 
is consistent with the real spectrum emission 
line correction; the dispersion of these distri- 
butions are around 0.05 dex. Therefore, we 
apply a systematic uncertainty of 0.05 dex 
to photometry corrected using the composite 
spectra. 

After emission line correction, the median 
photometric slope agrees well with the median 
spectroscopic slope in the UV band, with a 
value of —0.43. In the upper panel of Fig. [51 
we plot these two slopes for objects with SED 
fits. Most data points are distributed along 
the line representing ai,^ph — a,y^sp with con- 
siderable scatter. We believe that this scat- 
ter is mostly caused by low redshift quasars 
(typically z < 0.8) in which the 2200 A "hne- 
free" rest-frame continuum point is not cov- 
ered by SDSS spectra. In these cases, the un- 
certainty of the UV spectroscopic slope a^.sp 
may be larger than for high redshift quasar 
spectra. The value of a,y,sp can be more ac- 
curately measured for objects at higher red- 
shifts, which explains why the dispersion of 
Q^i/.ph — CKjy.sp is much smaller (bottom panel, 
Fig. [BJ . In the description of correlation and 
regression analysis, we consistently use auv 
to represent a^.ph because the data used to 
fit the slope is taken simultaneously with the 
X-ray data. 



10 



3.4. Photometric Shift 

The photometric data in SDSS and 2MASS 
are not observed simultaneously with the 
Swift UVOT data, so quasar variability can 
make flux comparisons uncertain. To make 
best use of these data, we must create some 
pseurfo-simultaneous data points by shifting 
the observed flux levels to be consistent with 
the Swift UVOT photometry, assuming the 
UV/optical SED shapes remain unchanged. 
With additional data points, we can place 
tighter constraints to UV slopes and lumi- 
nosities while still, in some sense, maintain- 
ing the simultaneous property of the dataset. 
The photometric shift follows strategies de- 
scribed below (see Table [5] and Fig. [7]). 

1. Because all quasars have flve simulta- 
neous SDSS photometric measurements 
we always interpolate or extrapolate the 
SDSS photometry to obtain fluxes at 
the UVOT fllter effective wavelengths, 
then shift aU SDSS and 2MASS pho- 
tometry to match UVOT. Because there 
is no overlap between the 2MASS and 
UVOT bands, it is impossible to shift 
2MASS photometry separately so we 
use the same amount of shift as the 
SDSS data. The 2MASS photometry is 
only used for a few high redshift quasars 
so the time difference between 2MASS 
and SDSS data wfll not affect the UV 
spectral slope for a majority of objects 
in our sample. 

2. We prefer to match in the UVOT U 
band. The Swift U band is very close to 
the SDSS u band, which signiflcantly re- 
duces uncertainties introduced by inter- 
polation or extrapolation. The U band, 
being bluer than the V and B bands, is 
also less contaminated by host galaxy 
light. 

3. If Swift U band photometry is unavail- 
able, we interpolate the SDSS u and g 
bands to match the Swift B band. If 



neither the U nor the B band is avail- 
able, we interpolate the SDSS r and g 
band to match the Swift V band. 

4. If none of U, B, or V bands is available, 
we extrapolate the SDSS g and u band 
to match one of the UVWl, UVM2 
or UVW2 bands, using the available 
UVOT filter with the longest effective 
wavelength. Matching these filters is 
only done as a last resort; the significant 
extrapolation required inevitably in- 
troduces considerable photometric un- 
certainties. Furthermore, these three 
bands frequently lie in the Lyman for- 
est where the quasar SED suffers from 
severe intrinsic and intervening absorp- 
tion (e.g.,|Rauch 19981). The continuum 
in this spectral region cannot be ap- 
proximated as a single power-law. 

5. When selecting the matching filters, we 
require the SDSS filters and the UVOT 
filters to both fall in the UV range (be- 
tween Lya and 5600 A) or within the 
EUV range, simply because the UV 
power-law cannot be extended to EUV 
region. 

Examples of photometric shifting results are 
presented in Fig. El 

3.5. Error Analyses 

For UVOT photometry, we adopt the 
photometric uncertainties produced by uvot- 
source and perform error propagation assum- 
ing these errors follow a Gaussian distribu- 
tion. The uncertainties of the SDSS and 
2MASS fluxes are calculated based on mag- 
nitude uncertainties in the SDSS DR5 quasar 
catalog. The frequency bandwidths plotted 
on the SED are converted from the corre- 
sponding FWHM of each filter. 

For the XRT data, we use the parameter 
uncertainties produced by the XSPEC error 
command. Because the X-ray photon counts 
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follow a Poisson distribution, which is unsym- 
metric, we calculate the upper and lower error 
bars separately. 

3.6. SED Models 

Because of the deficiency of the EUV 
data, the exact quasar bolometric luminosi- 
ties strongly depends on the model used to 
fit the BBB. Traditionally, emission in this 
"gap" is represented by a po wer-law contin- 
uum with a slope aox (e.g., iRichards et al 
20061 ). In this work, we consider a model 
which provides an upper limit to the flux of 
the BBB in the canonical case. This shape 
is inspired from the presence of a soft X-ray 
exc ess over a flat X-ray component reported 
by lArnaud et aP (|l985h in EXOSAT spec- 
tra of Seyfert 1 galaxy Mkn 8 4 1 and fur- 
ther studied bv I Walter fc FinkI (|l993l ) and 
Gierlinski fc Pond (|2004f ). The latter study 
found that this soft X-ray excess can be well 
fit by a black body of energy 0.1 - 0.2 keV. 
This result motivates us to constrain this fea- 
ture using a "bump" shape. In the following 
context, this model will be called the expo- 
nential decay model (hereafter EXP model). 
In addition, we use another model that di- 
rectly connects the high energy limit of the 
UV power-law and the lower energy limit of 
the X-ray power-law at 0.3 keV. We call it 
the triple power-law model (hereafter TPL 
model). In either model, total flux is the 
result of a UV and an X-ray component. 



Ua 



fu.VV + /i/,X- 



(8) 



EXP model In this model, the two com- 
ponents are 

/,,uv = 10^"-i.""-e('''''/'=^-)'"" (9) 

/^_x = 10^^:/"^e(-''''/'=^^)""' (10) 

The UV component is a power-law multiplied 
by an exponential decay term. The power-law 
slope Qfuv and scale factor /3uv are obtained 
by fltting photometric data points covered in 



the UV region (from 5600 A to Lyg). Canon- 
ically auv ~ —0.4 (e.g., IVanden Berk et al 



20011 ). The combination of exponential and 



power-law terms create a "bump" in the EUV 
region. The value of 7uv controls how quickly 
the flux decays given a value of Tb which is 
obtained from the SED fltting. We flx 7uv 
to be 1.5 so that at j/ < v(Lya) the contribu- 
tion from the exponential term is negligible 
and the total curve agrees well with the UV 
power-law given the best fit of T g. For exam- 
ple, i f we adopt 7uv = 1 (e.g., iGrupe et al 
20101 ) , a typical value of Tb '^ 4 Ryd from the 



SED fitting produces a discrepancy between 
the UV power-law and the total SED curve 
by ^ 0.1 dex at i^{Lya). 

The X-ray component is also a power-law 
multiplied by an exponential decay term. The 
power-law slope ax and scale factor /3x are 
obtained by fltting the X-ray energy spectra 
using XSPEC. The X-ray decay energy Tx 
is fixed at 0.3 keV. We adopt 7x = -8/3 
which assume s neutral hydrogen absorption 
(|LongaiiJll9"92| ). 

We emphasize that both of the UV and 
X-ray components are only mathematical 
expressions. Specifically, the value of Tb 
does not reflect the accretion disk temper- 
ature. Our goal is to place a reasonable up- 
per limit to the strength of the BBB which 
contributes a signiflcant fraction of the bolo- 
metr ic luminosity ( e.g. , iMathews fc Ferland 



1987: Zheng et al. 1997 



Teller et al. 20021: iShang et al 



Laor et al.l Il997 
2005|). There 



fore, the speciflc mathematical form is not 
important. 

Because the XRT observes photons be- 
tween 0.3 and 10 keV, these measurements are 
shifted to a higher rest energy range for high 
redshift quasars, i.e., we lack the soft X-ray 
data points for these objects. This affects our 
SED fltting because if we only flt the observed 
data points, the decay energy of the bump 
will be shifted to a higher energy band as the 
redshift becomes higher, which leads to addi- 
tional flux contribution. To solve this prob- 
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lem, we supplement artificial data points by 
extending the X-ray power-law from the mini- 
mum energy of the real data down to 0.3 keV. 
We use 0.3 keV as the lowest energy for the X- 
ray power-law based on previous X-ray stud- 
ies of quasars. These investigations found 
that the X-ray power-law c an extend from 
0.1 keV to > 10 keV 
19891: 



Laor et al. 



5-^ 



Turner fc Pound; 



200 



.. .19971: IGeorge et al. 

Page et al.ll2005l: lYoung et al.ll2009l ). We then 
fit all the data points with this SED model 
using the Levenberg-Marquardt algorithm, 
leaving Tb the only free parameter. The me- 
dian value of Tb is 5.96 Ryd. 

TPL model In this model, the UV and 
X-ray power-laws are fit in the same way as 
the EXP model. Instead of multiplying each 
power-law with an exponential decay term, 
we simply connect the UV power-law at Lya 
and the X-ray power-law at 0.3 keV and de- 
note the slope as auvx- The distribution of 
Q^uvx is presented in Fig. [S) 

In a canonical situation, i.e., auv= —0.4 
and ax= 2, the EXP model generates a bump 
in the EUV region which produces an upper 
limit while the TPL model yields a lower limit 
of the strength of the BBB. This simple anal- 
ysis may not apply to strong X-ray absorp- 
tion quasars. An example is shown in the 
first SED plot in Fig. [101 In this case, the 
absorption leads to a flat UV slope in addi- 
tion to low X-ray emission in soft X-ray band. 
In these cases, the EXP model does not pro- 
duce a "bump" but a "dip" in the EUV region 
and the corresponding integrated bolometric 
luminosity is not reliable. These quasars are 
flagged as "red" and will not be included into 
our cleaned sample defined below. Examples 
of fitted SEDs using these two models are pre- 
sented in Fig. [ini 

3.7. Bolometric Luminosity and Black 
Hole Mass 

To obtain black hole masses, we measure 
broad emission line widths with SDSS spec- 



tra. We employ the same software package 
used in previous SDSS spectral analyses (e.g., 
Wu et al. 2009; Vanden Berk et al. 2011, 
in preparation). The Galactic reddening cor- 
rections to all the spectr a are performed us - 
ing the extinction curve of lFitzpatrickl (jl999l ). 
Values oi E{B — V) are calculated following 
Schlegel et al.l (jl998l ). We use three compo- 
nents to fit a spectrum: a single power-law 
in the UV band (between the Lya emission 
line and ^ 5500 A), th e small blue bump, the 
UV ir on template from lVestergaard fc Wilked 
( 2001) and the optical iron template from 



Veron-Cettv et al.l (|2004l ). Each broad emis- 
sion line is fit by a single or multiple Gaussian 
profiles. To ensure that our software pack- 
age produces consistent results with previous 
studies, we compare our bla ck hole masses to 
these in lShenet aD(|2008bl). 

Following IShen et aD (|2008bl ). we adopt 
different black hole estimators depending on 
quasar redshifts. At z < 0.7, we use II/3 
plus Za(5100A), at 0.7 < z < 1.9, we use 
Mg II plus Za(3000 A), and at z > 1.9, we 
use C IV plus ;a(1350 A). The black hole 
mass is calculated using the following equa- 



tion (jShemmer et al.ll2006Hshen et al.ll2008bl) 



log Mbh = a + blog (ALa) + 2 log FWHM 

(11) 
in which A^bh is in solar mass M©, XL\ is 
in lO''^ ergs s~^, and FWHM is in km s~^. 
The two coefficients, a and b, in this equa- 
tion, are (a, 6) — (0.66,0.53) for quasars with 
z < 0.7, (a, 6) = (0.505,0.62) for quasars 
with 0.7 < z < 1.9, and (a, 6) = (0.672,0.61) 

for quasars with z > 1.9 (see Table El 

Mc Lure fc JarvisI l20o3 iMcLure fc DunlopI 



120041 : IVestergaard fc Petersonll2006( ^ 

For comparison purposes, we calculate the 
bolometric luminosity using the wavelength- 
depend ent bolometric correc tion (BC) fac- 
tors by iRichards et al.l (|200d ) which are ob- 
tained using the comp osite SED for a sample 
of SDSS DR3 quasars (JAbazaiian et al.ll2005 : 
Schneider et al.l 120051) , so that iboi sdss = 
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BC • Xl\. These are bolometric luminosities 
without using simultaneous observations and 
does not include the BBB component, which 
differs from the bolometric luminosity we ob- 
tained by SED integration. 

We then compare some importan t cor- 



relation s and distributions with Slien et al 
( 2008b[ ). The black hole masses vs. bolo- 



metric luminosities of 923 quasars are plotted 
in Fig. Iljn . The locations of these quasars 
are consisten t with the location s of quasars 
in Fig. 11 of IShen et al] (|2008bl) . indicating 
that most quasars are accreting at a sub- 
Eddington level, i.e., 0.1 <iboi/-^Edd< 1- 
The few quasars that fall above -/jboi/-^Edd= 1 
(they are accreting at super-Eddington level) 
are narrow lin e Seyfe r t 1 (NLSl) galax- 
ies (|Shen et al.l l20Q8al: iGrupe et all |2010[) . 
The median values of logFWHM/(km s~ ) 
are 3.62, 3.63 and 3.75 for H/?, Mg ii and 
C I V, respectively, w hich are consistent 
with IShen et al.l (|2008bl ). By comparing the 
FWHM calculated using both of H/3 and Mg ii 
or Mg II and C iv we find a potential source 
of bias in the calculation of black hole masses 
as the C iv FWHMs are systematically larger 
than those of H/3 and Mg ii. However, we still 
use the black hole masses estimated by C iv 
because this offset is small compared to the 
dispersion of the black hole mass distribution. 

We also calculate the bolometric luminosi- 
ties by integrating the two SED models from 
5600 A to 20 keV. In general, emission in this 
region contributes the majority of the quasar 
luminosity for typical Type 1 quasars. 

In Fig. [121 we show the differences be- 
tween the bolometric luminosities measure- 
ments by integrating the two SED models, 
and the bolometric luminosities calculated 
using BC. The EXP model and the TPL 
model produce bolometric luminosities av- 
erage higher and lower than the iboi from 
the BC, respectively. Because the bump 



^The others are not displayed because of low SDSS 
spectral quality. 



emission from the EXP model is more sen- 
sitive to the UV and X-ray observed spec- 
tral shapes than the TPL model, it has a 
larger dispersion. A small fraction of the 
EXP models produces less bolometric lu- 
minosity (logLbol,EXP - logibol,SDSS < 0) 
than the TPL models. Some of these ob- 
jects are quasars with strong intrinsic ab- 
sorption (the black shaded area in Fig. [T2|) . 
The reason why the TPL model produces less 
bolometric luminosity on average than the 
BC meth od is because the bo lometric cor- 
rection bv [Richards et al.l (|2006l ) includes the 
infrared wavebands, which contribute about 
37% of the total flux. After correcting for 
contributions from this component, the en- 
tire histograms of both the EXP and the TPL 
models should move positively by ~ 0.2 dex. 
Consequentially, the TPL model on average 
produces a bolometric luminosity consistent 
with that from the BC and the EXP model 
on average over produces iboi by ^ 0.3 dex. 

4. Catalog Description 

Our final catalog contains 843 quasars with 
0.0129 < z < 4.5766 and -30.24 < Mi < 
—22.01. There are 675 objects observed by 
UVOT, 805 observed by XRT and 637 ob- 
served by both UVOT and XRT. The X-ray 
detection rate of the entire catalog is sa 60%. 
Among objects in the catalog, 345 objects are 
detected by both UVOT and XRT, so that we 
are able to determine their SEDs from UV 
to X-ray bands. These objects constitute our 
parent sample to evaluate the flux contribu- 
tion from the BBB. This catalog contains pa- 
rameters directly measured from UVOT and 
XRT data and quantities derived based on 
these measurements such as black hole masses 
and bolometric luminosity. Columns in this 
catalog are described in Table |7l Comments 
for special individual objects in this catalog 
are listed in Appendix [71 Because the cata- 
log contains 63 columns, we publish the full 
catalog in electronic format. In Table |8l we 
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represent ten columns of some objects in our 
catalog. 

5. Correlation Analyses 

In this section, we investigate the aox~ 
L ' relationship using a selected sam- 

2500 A f o 

pie of Type 1 quasars from our catalog. 
This correlation has been described in a 
se ries of paper s ( Ayni fc TananbaumI 1982 ; 



Marshall et al.' 1984; Tananbaum et al. 



Anderson fc Margon.l987.:.Wilkes et al 



2003HStrateva et al.ll200a ISteffen et aLll2006t 



1986; 



1994; 



Pickering et al.ll994lAvni et al.ll 99511 Vignali et al 



Just et ahlbOOTl) ; as quasars become more lu- 
minous, their SEDs from UV to X-ray bands 
become softer, i.e., less X-ray emission with 
respect to UV emission. In the study of 372 
type 1 quasars in Just et al. (2007; hereafter 
J07) , the UV/optical data ar e drawn from the 



SPSS DR3 quasar catalog ([Schneider et al 
200l, lUE (JBoggess et aIlll978l : lKondo et al. 



1989^ and the COMBO- 17 survey (^Wol f et al 
20031 which covers the E-CDF-S (Wolf et al 
20oS i. The X-ray data are from ROSAT, 
Chandra, and XMM-Newton observations. 
Because the time between UV/optical and 
X-ray observations could span a time scale 
of years, quasar variability inevitably intro- 
duces scatters to the a„x-L = relation. In 

2500 A 

this section, we investigate whether this scat- 
ter can be reduced by simultaneous UV and 
X-ray observations. 

5.1. Sample Selection 

In this work, we use 637 quasars with si- 
multaneous observations by both UVOT and 
XRT. Because this sample has a relatively low 
X-ray detection rate ('^ 65%, see below) we 
applied an X-ray exposure cut-off to the XRT 
observations. The goal is to obtain a homoge- 
neous, optically-selected sample with a rela- 
tively high X-ray detection rate while retain- 
ing a sufficiently large sample size. The sam- 
ple size and X-ray detection rate as a func- 
tion of exposure cut-off is presented in Fig.lTSl 



The X-ray detection rate is ^ 65% if we use 
the entire parent sample. These censored 
data can be han dled by the A SURV software 
package (Lavall ev et al.lll992l l. but the large 
fraction of non-detections makes the results 
unreliable. On the other hand, if we increase 
the exposure cut-off, the detection rate can 
rise to > 95%, but only a few objects remain 
in the sample and the statistics are very poor. 
In this paper, we use two samples to study the 
aox~L ' relation: the large catalog sam- 
pie which contains all qualified quasars se- 
lected from the parent sample with no XRT 
exposure cut-off; and the small catalog sample 
which is derived from the large catalog sample 
except that we apply an XRT exposure cut- 
off of 10 ks. The large sample has a sample 
size of over 400 quasars; the small sample has 
a higher X-ray detection rate (~ 85%, see Ta- 
ble [T0|). In the following analysis, we supple- 
ment both o f these samples with AGNs from 
Grupe et al.l ( 20101 ). which contains data for 



88 AGNs observed simultaneously by Swifi^. 
We define this sample supplemental sample or 
the GIO sample. The GIO sample is soft X- 
ray selected and is composed of low luminos- 
ity Type 1 AGNs at low redshifts {z < 0.4). 
Because the GIO sample is selected in a differ- 
ent way from our sample, this sample is not 
merged with our samples, but is only used for 
comparison purposes. 

The catalog parent sample is mixed with 
different types of quasars. To obtain a clean 
catalog sample, we exclude objects that fall 
into any of the four categories listed below. 
The quasar types are determined by explor- 
ing previously published studies and by cal- 
culating relevant quantities, e.g., radio loud- 
ness. The numbers of quasars excluded in 
each type are tabulated in Table [H Because 
the GIO sample was not selected in the same 
way as our catalog sample, we only use them 



^ There are 92 objects in total in the iGrupe et al.l 
1 I2OIOI ) study, but 4 objects were not observed with 
UVOT photometry. 
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to represent simultaneously observed objects 
at relatively low redshift and low luminosity. 
Therefore, the exposure time cut-off is not ap- 
plied to the GIO sample. The GIO sample is 
X-ray selected, so all objects are detected by 
XRT. 

RL quasars They frequently have sub- 
stantial X-ray flux from the radio jet, which 
leads to a higher aox than thos e of radio quiet 



quasa rs ( Worrall et al. 1987 : Brinkmann et al. 
2000h . 



We adopt the "radio loudness" (R*) de- 
fined by the ratio of monochromatic lumi- 
nosities at rest-frame 5 GHz and 2 500 A 
(IStocke et"all Il992l: Idella Ceca et al. 



Lu et al.ll2007t iMiUer et al.ll2009l) . e.s.. R* 



19941: 



L 



GHz/-^, 



5 L;J1z/""2500 A' 

For most of the quasars in the parent 
sample, we calculate the fc-corrected Lg GHz 
from the 20 cm peak flux from t he Very Large 



Array (VLA) FIRST S urvey (jBecker et al. 



199alWhiteet al.lll997l) listed in the SDSS 
assuming an average 



DR7 quasar catalog, 

value of radio spectral i ndex ar 



Kelle rmann et al.lll989t iKomossa et al.ll2006t 
Luet al. 200^. 

For objects not detected by the FIRST 
survey (but still covered), we estimate ra- 
dio luminosity upper limits using the sensi- 
tivity limit of the FIRST survey, which is 
1 mjy. This provides upper limits of R* . If 
logi?* < 1, we classify this object as RQ. If 
log i?* > 1 , we check other resources of radio 
surveys, such as the NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998 ). 

For objects not covered by the FIRST sur- 
vey, we examine 15' NVSS images for strong 
nearby radio sources. If we cannot confirm 
the quasar is RQ, we conservatively exclude 
it from the clean sample. 

We preferentially adopt L ° from SDSS 

spectral fitting results, unless 2500 A is not 

cal- 



covered, in which situation we use L 
culatcd from the SED fitting. 
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Lensed quasar The fluxes from these 
quasars are amplified, so these obj ects are re- 
move d from the clean sample (e.g.. ljust et al. 



20071 ). 



Blazars Blazars, including BL Lac objects 
and FSRQs, usually have large amplitude 
varia tions in UV/optical (e.g., I Rani et al.l 
20101 ) as well as strong radio and X-ray emis- 



sion f e.g. jBeckmann et al.ll2003l : lPadovani et al 
20071 ). As a result, the values of Ofox can 
vary significantly over time. The UV/optical 
spectra of BL Lac objects usually have a 
featureless continuum without any emission 
lines (e.g.. [Blandford fc Reeslll978l: iKollgaardI 
119941 : IPlotkin et al.l 120081 : lAbdo et al.ll2010[ ). 
Many of them are RL and their emissions are 
believed to be relativistically b eamed from 
the jet (jUrrv fc Padovanil Il995l ). Although 
some FSRQs exhibit similar broad emission 
line features as Type 1 quasars, their radio 
emission is still beamed and variable. All ob- 
jects classified as blazars are also excluded 
from the clean sample. 

Reddened quasars Our sample contains 
a number of quasars with shallow UV/optical 
slopes and/or strong soft X-ray decline. The 
relatively fiat UV/optical s pectra are at 



tributed to du s t red d ening (IRichards 



20011: 



20041 ) 



[Richards et"al1 120031 : iHopkins et al 
although observations of individual objects 
suggest that son ie slopes could be intrinsi- 
cally steep (e.g.. iHall et al.l 120061 ). Because 
the gas density of the AGN BLR is much 
higher and radiative transfer effects are not 
the same as in the low density regime (e.g., 
iHao et al.l 120051 ) . traditional approaches to 
correct reddening by the Balmer decrement 
cannot be applied for quasars. In this work, 
we only flag these objects and exclude them 
from the cleaned sample. 



Following [Richards et al.l (|2003l ). we use 



the relative color to define the "dust-reddened" 
quasars. The relative color is calculated by 
subtracting the median colors of quasars at 
the redshift of each q uasar from its m e asure d 
colors. As argued bv iRichards et al.l (|2003[ ). 
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the quantity A (g — i) = {g—i)—Mcdia,ii{g^i)z 
is an excellent redshift-independent surrogate 
for the photometric spectral index. In Fig.[T4j 
we plot A{g — i) vs. redshift for all the objects 
in our raw catalog. Objects to the right of the 
dashed line are flagged as "dust-reddened" 
quasars. These quasars, which comprise 8.5% 
of the raw catalog, are excluded from our 
clean sample. 

Broad absorption line (BAL) quasars 

BAL quasars are excluded from the clean 
sample because the emission line absorption 
is found to be associated with the contin- 
uum absorption in both of the UV/optical 
and the X-ray band s fe.g-. lMathur et aUlQQSl : 



Brandt et al.l 120001 ). This property can sig- 



nificantly affect the spectral indices and flux 
values obtained by fltting photometric data 
points, which produces inaccurate values of 
Qox- Most BAL quasars in the parent sample 
are identified in th e BA L quasar ca talogs by 
Trump et al] (|2006l) and lGiJjson eTa l. (2009). 
We found 50 BAL quasars in the parent sam- 
ple; seven of them are also identified to be red 
quasars. 



5.2. The ttox L 



A 



Relation 



In this section, we compare the dispersions 
of the aox-L ° relation between the J07 

2500 A 

sample and our cleaned catalog sample. We 
show that the dispersion can be reduced using 
simultaneous UV and X-ray observations. Be- 
cause the exact values of the dispersion could 
be method-dependent, we use two methods 
in our analyses. We will present two sets of 
results based on both of the small and large 
cleaned catalog sample as well as for the com- 
bined sample. 

Because our sample contains censored 
data, we use the ASURV software pack- 
age (JLavallev et al.l Il992.) to perform sta- 



tistical analysis. This package includes the 
Expectation-Maximization (EM; Dempster et 
al. 1977) and the Buckley- James (BJ; Buckley 
& James 1979) algorithms, which we can use 



to perform linear regression and dispersion es- 
timation. The methods differ in that the EM 
algorithm estimates the residual assuming a 
Gaussian distribution while the BJ algorithm 
assumes the Kaplan-Meier distribution. The 
correlation and regression results are tabu- 
lated in Table [11] We discuss the results and 
their implications below. 

The clean catalog and the combined sam- 
ples both exhibit strong correlations between 



ttox and L 



although the absolute val- 
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ues of the correlation coefficients are slightly 
lower than the J07 sample. The relatively 
large fraction of undetected objects in the 
large clean catalog sample smears the correla- 
tion, but it is more evident in the small clean 
catalog sample with a higher X-ray detection 
rate. 

The slopes (— 0.16±0.02 for the small sam- 
ple and —0.15 ± 0.01 for the large sample) 
are both steeper than the J07 sample. The 
intercepts are also larger, but they are both 
consistent within 2(t uncertainty. The com- 
bined sample has a shallower slope than the 
clean catalog sample; this change is caused 
by the GIO sample. The possibility that 

the a^r-L ° relation is non-linear is pro- 
2500 A ^ 



pos ed in the study of l^yilkes et al.l (|l994l ) 
and [Anderson et al.l (12003). As discussed by 



Wilkes et al.l (199J), the difference in slopes 



is likely caused by the varying host galaxy 
contribution to the L ° measurement at 

2500 A 

low redshift. This may be the reason for 
the shallower slope of the GIO sample since 
no host g alaxy contribut io n cor rection is ap- 



plied Q. IStrateva et"all (|2005l) found that 



their sample does not offer significant evi- 
dence for a non-linear a^v-L ° relation. 

™ 2500 A 

Although they obtained a shallower slope of 
— 0.09±0.02 for the low luminosity (log Lyv < 
30.5) and -0.13 ± 0.02 for the high luminos- 



^We do not correct for host galaxy light in our sample 
either, but since the majority of our objects are lu- 
minous quasars at higher redshift than GIO, the host 
galaxy contamination is much smaller. 
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ity (logLuv > 30.5) sample, they argue that 
the difference in slopes is likely an artifact of 
the addition of five outlier AGNs, which are 
probably X-ray-absorbed Seyfert galaxies at 
z < 0.22. From Fig. [TH we do not see sus- 
pected low luminosity "outliers" , but we still 
cannot exclude the possibility that the shal- 
lower slope is caused by host galaxy contam- 
ination. In addition, we also note that the 
J07 sample contains a number of the most 
luminous quasars with logLuv > 31.5, but 
excluding these high luminosity quasars does 
not reduce the slope significantly (see regres- 
sion results of J07t[1 in Table [nj). A careful 
removal of host galaxy contribution to the im- 
ages is probably nece ssary for further checks , 
e.g., by using GalFit (|Peng et al.ll200ll2010l) . 
From Table [11] it is clear that the stan- 
dard deviation of the clean catalog sample is 
smaller than the J07 sample by about 13%- 
19%. The large clean sample has the largest 
dispersion (BJ algorithm), which is compara- 
ble with the J07 sample, but this large dis- 
persion is mostly caused by a large fraction 
of undetected X-ray sources. Using the EM 
algorithm, the J07 sample still has the largest 
dispersion among all samples. The small cata- 
log sample exhibits a dispersion reduced even 
more by 18%-25%, compared with J07. In 
general, the combined sample has an even 
smaller intrinsic dispersion because the GIO 
sample has a higher degree of simultaneity. 

The dispersion of the a„-x-L ° relation 

. . . . 2500 A 

does not exhibit a luminosity dependence (see 
Figure 7 of Just etal. 2007). From Fig. [T5j we 
do not see the dispersion showing any evident 
dependence on luminosity (at least for X-ray 
detected quasars), either. 

5.3. Correlations Between UV/optical 
and X-ray Spectral Indices 

Previous studies found that AGNs with 
bluer optical/UV spectra have relatively 



steeper X-ray spec tra ([Walter fc FinkI 1 1993 : 
Grupe et al.l l2008l) . U sing simultaneously 
observed nearby AGNs, iGrupe et al.l (|2010l ) 



found a mild correlation between auv and ax 
(See Figure 10 of their paper). In Fig. [TCI we 
plot a\jy from UVOT photometric data vs. 
the X-ray spectral slope a^ for objects from 
our clean catalog sample (objects whose a^ 
were fixed to —1 are excluded). The Spear- 
man correlation coefficient for the data is 
Ps — — 0.136 (Po — 0.125), which indicates a 
very weak correlation but the value of Pq in- 
dicates low confide nce leveO- As claimed by 
Grupe et al.l ( 2010l) . the correlation between 
auv and a^ is primarily driven by BLSl with 
ax> —1.6. By excluding NLSl and ax< —1.6 
spectra objects from our sample, we find that 
the correlation coefficient Ps = —0.22 with 
Pq = 0.016, which still does not indicate a 
significant correlation. 

To investigate whether the lack of this 
correlation in our sample is due to the 
wider range of re dshifts compared to the 



Grupe et al.[ (|2010() sample, we select the low- 



redshift counterparts from our sample with 
z < 0.4, marked in blue in Fig. [THI These 
objects exhibit a stronger correlation, with 
Ps = -0.468 (Po = 0.058). Although this 
subsample contains only 17 objects, we ar- 
gue that this result is expected and due to 
changes of measured spectral slope in soft and 
hard X-ray bands. In the GIO sample, most 
AGNs are at low redshift so the rest-frame 
X-ray spectra cover both soft and hard X-ray 
bands. Because this sample is soft X-ray se- 
lected, objects in this sample usually exhibit 
strong soft X-ray emission with respect to 
hard X-ray emission, which is caused by the 
soft X-ray excess. The majority of our sample 
is composed of quasars with much higher red- 
shifts, so the rest-frame X-ray spectra cover 
relatively less soft X-ray band. As a result. 



°The low luminosity sample logL 



< 31 in J07. 



^The null hypothesis here is that the correlation does 
not exist, so the lower Pq is, the more confident we 
feel on this correlation and vice versa. 
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the measured photon indices are mostly based 
on hard X-ray data. 

We further investigate whether the corre- 
lation differences are due to the differences 
between the luminosity range in our sample 
and that of GIO. As seen in Fig. [T5l the 
GIO sample is composed mostly of low lumi- 
nosity AGNs with a luminosity upper limit 
of loeL i ~ 30.5. We select a subsam- 

° 2500 A 

pie of quasars from the clean catalog sample 

with logL ! < 30.5 (green diamonds in 

° 2500 A ^^ 

Fig. [T6| and found that the correlation co- 
efficient is only 0.043 with a very low confi- 
dence level Pq — 0.726. Therefore, our data 
do not support the argument that the differ- 
ences between axi v-Qy correlation found by 
iGrupe et al.1 (|2010l ) and in our work are due to 
luminosity differences, but more objects are 
needed to verify this. 



5.4. 



Correlation between ciox and Spec- 
tral indices 



Usi ng the low redshift AGN sample. iGrupe et al 
( 2010r i found a correlation between aox 
and Qix, i.e., AGNs with softer X-ray spec- 
tra tend to be X-ray weak relative to UV 
band, which is consiste nt with the results of 
Atlee fc Matlnul (|2009l ). In Fig. [17l we plot 
Qox vs. ax for 129 objects in our clean catalog 
sample. We do not see a significant correla- 
tion with ps = -0.062 (Po = 0.484). Similar 
to the method described in Section 15. 3[ we 
select a subsample with z < 0.4, which is the 
low redshift counter part of the GIG sample. 
The Q!ox and a^ of this subsample exhibit a 
strong correlation, although this subsample 
contains only 17 objects. 

To see if the lack of correlation is due to 
redshift, we plot the high redshift (z > 1.5) 
objects of the catalog sample in red. These 
objects are at a different location from the 
GIO sample with relatively higher values of 
Qx- As we argued in Section 15.31 the mea- 
sured value of spectral slope may depend on 
the rest-frame energy range. Because the ob- 



served energy range is fixed, we do not see 
the soft X-ray energy spectra for high redshift 
quasars and thus photon indices measured for 
these objects suffer less from the soft X-ray 
excess. 

We also display the low luminosity sub- 
sample (logL ! < 30.5) using diamonds 
in Fig. [T71 which contains only eight objects. 
The subsample size is too small and we can- 
not decide if lack of correlation is related to 
the observed luminosity range. 

Compa red with th e posit ive correlation 



found in iGrupe et al.l (J2010[ ), we plot the 



margi nal linear correlation found bv lYoung et al 
()2009l ) in Fig. [T71 which show s a negative 
correl ation. The regression by lYoung et al. 



(|2009l ) in gener al agrees with our d ata trend. 
The sample in I Young et al.l ( 20091 ) contains 
RQ quasars over a redshift range of z = 0.11- 
5.41, which is very similar to our s ample . 
However, both the lAtlee fc Mathuij (|2009t ) 
and the GIO sample contain soft X-ray se- 
lected AGNs at low redshift. Because the 
fixed observed energy range, the spectral in- 
dices we measured for high redshift quasars 
represent the hard X-ray spectral shapes and 
are less affected by the soft X-ray excess which 
produce lower values of ctx ■ The marginal 



correlation found in lYoung et al.l (20091) is 



not seen in our sample is likely because they 
fit their spectral over the 0.5-10 keV band 
which covers less portion of the soft energy 
band and are less vulnerable to the soft X- 
ra y excess. Therefore , the correlation found 
by JGrupe et al.l ()2010l ) basically implies that 
AGNs with softer SEDs over the UV and X- 
ray bands tend to have stronger soft X-ray 
excess. This is consistent with the argument 
that the soft X-ray excess is a tail of the 
BBB in the EUV band, which originates from 
the thermal emission from the accretion disk. 
The lack of correlation in our sample and 
margi nal correlation found by JYoung et al 



( 20091 ) implies that the hard X-ray genera- 
tion process is relatively independent of the 
process producing the BBB photons. 
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Grupe et al.l (|2010l) also found a correla- 
tion between aox and any, which is not seen 
in our sample. The Spearman correlation co- 
efficient is ps = 0.166 (Fo = 0.014). The low 
redshift subsample {z < 0.4) contains only 
nine objects which reside at similar location 
as the GIO sample and appear to exhibit a 
strong correlation with ps = —0.683 (Pq = 
0.042). However, the objects in our sam- 
ple with higher redshifts cover a much wider 
range of Oox values. 

The GIO sample is X-ray selected, so this 
sample will naturally include objects with 
stronger X-ray with respect to UV emission. 
Our sample is optically selected so it may con- 
tain objects with relatively weak X-ray with 
respect to UV/optical emission. It is clear 
to see from Fig. [15] that our sample contains 
objects with lower values of ctox- If we com- 
bine our sample with the GIO sample, we see 
a mild correlation between aox and auVi but 
with a large dispersion for high auv objects. 

5.5. Correlation Between iboi/-^Edd vs. 
Spectral Slopes 

It has been reported that the Eddington 
ratio Lha]/L-p M is correlated w ith auVj ctx, 



Grupe et al.ll2010l ). Intuitively, 



and ttox (e.g., 
this reflects that the accretion state affects the 
AGN spectral shape. Using the clean catalog 
sample, we find a mild correlation between 
the Eddington ratio and the UV spectral in- 
dex (Fig. [T9|) . The Spearman correlation co- 
efficient is Ps = 0.35 (Po < 10"^). The linear 
regression result is 

logLboi/l^Edd - (1.194±0.302)auv+(0.102±0.094) 

(12) 
In Fig. 1191 we over plot the linear regres- 
sion results for BLSls (magenta dash dotted 
line), NLSls (blue dash dotted line) and both 
(black dotted line). The BLSl linear relation 
is clearly more consistent with our result be- 
cause most objects in our sample are Type 1 
quasars. We also distinguish objects in our 
sample at different redshifts but we do not see 



any systematic offset for quasars at different 
redshifts. However, because the Eddington 
ratios used in Fig. [Tulare calculated using the 
bolometric luminosity under the EXP model, 
the values of iboi/^Edd are not entirely inde- 
pendent of auv- Thus, an independent mea- 
surement of Pboi/iEdd is required to verify 
this correlation. 

Compared to auv, the value of bolomet- 
ric luminosity is much less dependent on a^ 
(Fig. EO]) . The Spearman correlation coeffi- 
cient is Ps — —0.36 (Pq < 10"'^) which indi- 
cates a mild correlation. We perform linear 
regression to all data points in our sample 

logLboi/PEdd = (-1.187±0.211)ax+(-1.740±0.175) 

(13) 
The slop e is generally consiste nt with that 
found bv lShemmeret al.l (l2008l) (-0.9 ± 0.3) 
and shallower tha n the slope reported by 
Grupe et al.l (l2010l) (-1.65±0.26). InFig.[20l 



we distinguish objects at different redshifts in 
different colors. Again, quasars at low red- 
shift {z < 0.4) have relatively lower values 
of ax, because of the soft X-ray excess. Af- 
ter excluding these objects from our clean 
sample, the correlation coefficient is ps — 
—0.30 (Po = 0.001), and the linear regression 
slope is —0.88 ± 0.61. The enhanced correla- 
tion for low redshift quasars indicates that the 
accretion state change causes the change of 
the accretion disk temperature. Specifically, 
higher accretion rate increases the disk tem- 
perature, which leads to a higher level of soft 
X-ray excess. 

Grupe et al.l (|2010l ) also found a strong cor- 
relation between Lboi/iEdd and aox- This 
correlation was not s een in the sample of 
IShemmer et"aLl (|2008l) . In Fig. [HI we see 
a very weak correlation between these two 
quantities with ps = -0.16 (Pq < 10"^). By 
distinguishing quasars in different redshifts 
in colors, we notice that the strong corre- 
lation between Pboi/^Edd and aox stronger 
for low redshift quasars. Because of meth- 
ods used in sample selection, the GIO sam- 
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pie consists of mostly X-ray bright AGNs, 
thus relatively higher value of aox- Because 
of this bias, many UV/optical bright and X- 
ray normal/faint quasars are excluded. After 
including these quasars, the correlation be- 
comes much less significant. Even by combin- 
ing GIO sample with our sample, the correla- 
tion remains weak with significant scatter. 

5.6. X-ray Slope versus Redshift and 
Luminosity 

We do not find that the X-ray pho- 
ton index F has any significant correla- 
tions with redshift, UV or X-ray luminosity 
(Fig. [2D). The Spearman correlation coeffi- 
cients are p^ = — 0.14(Po — 0.003) for red- 
shift, ps = -O.f 3 (Pq = 0.13) for i » , and 
Ps = -0.09 (Fo = 0.27) for 2.3 keV°°^dding 
the GIO sample does not make these corre- 
lations stronger. This result confirms previ- 
ous s tudies with smaller or comparable sam- 



ples (IPage et al 



Shemmer et al 



2004 iRisaliti fc Elvid 



2005 



Vignali et al.l 



Kelly et al.ll2007l : lYoung et al.ll2009[ ). 



2005 



2005; 



6. Conclusions 

We have compiled an optically selected 
quasar catalog with serendipitous and simul- 
taneous UV/optical and X-ray observations 
with the Swift] observatory. The catalog is 
generated by matching the ~ 3.5 year Swift 
pointings from November, 2004 to June, 2008 
to the SDSS DR5 quasar catalog. For each 
object, the sky images observed by either 
UVOT or XRT are carefully selected to ensure 
high image quality. We derive the composite 
UVOT photometry and XRT energy spectra 
by stacking all archival data to generate the 
deepest sky images. The resultant SEDs re- 
flect the time-averaged shape of quasar emis- 
sion with simultaneous observations at multi- 
wavebands. The catalog contains 843 objects. 
There are 637 objects (~ 76%) that have 
UVOT and XRT observations, 168 objects 
(~ 20%) that only have XRT data and 38 



objects that only have UVOT data (- 4%). 
Among all the 675 objects with X-ray cover- 
age, 460 (~ 60%) are detected, which rises 
to 85% amongst source with at least 10 ks of 
XRT exposure time. We construct SEDs for 
all objects with both XRT and UVOT data. 
In this work, we focus on 637 objects with 
both X-ray and UV observations. We sup- 
plement UVOT photometry with SDSS and 
2MASS data if available. All the photometric 
points are corrected for the effects of emis- 
sion lines, and fluxes from SDSS and 2MASS 
are shifted to match the flux levels of Swift 
UVOT data. 

We fit SEDs using the EXP and TPL mod- 
els, attempting to constrain the flux contri- 
bution from the BBB. In most cases, the 
EXP models create a bump in the EUV re- 
gion, producing an upper limit on the BBB 
emission, while the TPL model connects the 
SED points at 1216 A and 0.3 keV, produc- 
ing a lower limit on the BBB emission. After 
correcting for the contribution from IR emis- 
sion, the TPL model produces bolometric lu- 
minosities consistent with those estimated us- 
ing BCs from composite quasar SEDs, while 
the EXP model produces bolometric lumi- 
nosities on average 0.3 dex higher than the 
TPL model. 

We identify two clean samples (large and 
small) selected from our catalog, and sup- 
plem ent each sam p le wit h 88 nearby AGNs 
from iGrupe et all (|2010t ). We r e-visit the 
aox--^„^„„ .° relation presented bv lJust et al. 



2500 A 



(|2007n . We use the EM and Buckley- James 

methods to compare the intrinsic scatters of 

the a„y:-L ° relationship of our and the 

2500 A ^ 

J07 sample. These two methods consistently 

indicate that the dispersion based on our sam- 
ple is reduced compared to J07 by 13% to 19% 
using the cleaned catalog sample and 18% to 
25% using the combined sample. 

Firmly establishing the ao^-L ° rela- 
tion in AGNs is an important step toward 
understanding energy generation mechanisms 
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of AGNs. Our work has verified the corre- 
lation again, and has shown that the corre- 
lation is even tighter after reducing or elim- 
inating scatter due to variability. An ad- 
ditional source of scatter in the correlation 
could be due to an intrinsic a„^-L ■> re- 

2500 A 

lation, which may have a different slope for 
each q uasar. In Fig. 7 of lVasudevan fc Fabian 
( 20091 ) ■ they present simultaneously and non- 



simultaneously observed SEDs in the UV and 
X-ray bands. The dramatic difference be- 
tween two SEDs for three of these objects 
indicates the variation of aox with time for 
a given AGN. This variation then reflects the 
change of accretion state of the central engine 
and a different slope for an intrinsic and for a 
global relation would produce scatter around 
the global relation. In order to see the in- 



trinsic Qfox variation with L 



for an given 
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AGN/quasar, it is necessary to gather long- 
term simultanous observational data. In our 
fugure work, we will perform time-resolved 
UV/X-ray data anlayses to selected targets 
from our catalog which was observed with 
high cadence, which will determine the contri- 
bution to the scatter of the global aox-L 

. . . 2500 A 

relation from the intrinsic variations of X-ray 
with respective to UV emissions. 

We also investigate correlations between 
spectral shapes in different wavebands and 
-Lboi/-^Edd; and compare the results found in 
GIO. Our low redshift {z < 0.4) counterparts 
to the GIO sample verify significant correla- 
tions exist between auv and ax, Q^ox and ax, 
aox and auv, iboi/-^Edd and auv -^boi/iEdd 
and ax, iboi/^Edd and aox, which physically 
implies that the accretion status plays a fun- 
damental role in shaping the quasar SED be- 
tween UV and soft X-ray band. This supports 
the argument that the BBB is produced by 
the disk emission and the soft X-ray excess 
is a result of thermal emission from accretion 
disk. However, for high redshift quasars, the 
measurement of X-ray spectral slope covers 
less soft energy band and is less affected by 
the soft X-ray excess. As a result, the cor- 



relations between spectral shapes are much 
weaker after including high redshift quasars. 
This implies that the hard X-ray emission is 
relatively independent of the thermal emis- 
sion on the accretion disk. 
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7. Comments on Individual Objects 

The following sources possess special spec- 
troscopic or photometric features or have un- 
usual classifications. 

SDSSJ021702.66-082052.3. This is a flat 
spectrum radio quas ar (FSRQ) cataloged by 
The optical spectrum 
continuum with very 



mined. 



Massaro et al 



(l2009f ) 
exhibits a featureless 
weak Mg ii A2798 and H/3 emission lines and a 
few narrow lines such as [O ill] AA4960, 5008. 
SDS SJ07i 1 1 0. 70+ 31 1 200. 2. This 



IS an 



FSRQ (|Lhi fc Zhan e 20()2: lHealev et al.ll2007h 
SDSSJ074 625. 86+ 2 54902.1. T his 



IS a 



high redshift FSRQ (IMassaro et al. 
with z = 2.979 and 
covered by Swift 



20091) 



an 



MeV blaz a r dis- 

(ISambrunaet al.l l2006t 



Tueller et al] l2008t Ijollev et al.l l2009l ) and 



thus in the BAT-s elected AGN catalog by 
Winter et al.l ( 20091). It i s also observed by 



Suzaku (jWatanabe et al] |2009|) . The UV 



spectrum contains broad emission lines such 
as Lya, C iv A1549 and Cm A1908, which is 
unusual for a blazar. 

SPSS J 081 331 . 28+ 2 54 503. 0. This is a con- 



firmed lensed quasar (IWilliams et al.l l2008t 



Congdon et"alll2010l:llnada et al.ll2oTc 

SDSSJ083148. 87+042939.0. This is an 
FSRQ with peak radio flux of ^ 1 Jy. Its 
optical spectrum exhibits a featureless con- 
tinuum with an extremely weak Ha emission 
line. 

SDSSJO8374O.24+245423.I. This object 
is a blazar a nd cl assified as an FSRQ by 
Healev et all (|2007l ). The UV spectrum ex- 
hibits strong Mg 11 and C in] A1908 emission 
lines. 

SDSSJ090821. 01+045059.4. The Mg 11 
emission line at ^2798 A is barely detected, 
H/3 is extremely weak and the two nearby 
[O III] lines are quite prominent, indicating 
that the BLR is obscured and only the NLR 
is seen. This object has been identified as a 



strong radio source (e.g.. lGriffith et al.lll995[ ) 



but the exact classification is not yet deter- 



SDSSJ092703. 01+390220. 8. This is an 



FSRQ with a strong rad io jet (jLiu fc Zhang 



2002 : Healev et al.l 120071 ). The UV / optical 



spectrum also contains prominent Mg 11 A2798, 
H/3 and H 7 emission lines. 

SDSSJO94215.I2+ 090015. 8 This is quasar 
with double-peaked Ha. The FWHM of 
the broad component is extremely wide and 
reaches ^ 40, 600 km s''^, wh ich is the broad- 
est known (jWang et al.ll2005[ ). The broad Ha 
suggests that the emission re gion is close to 
the b lack hole r 
200i). 

SDSSJIOI4O5. 89+000620.3. This quasar 



80-100 r„ ( Stratevaet al. 



i s clas sified as a Seyfert 1.8 by iDong et al. 
(|2005l ). The optical spectrum exhibits 



strong and double-peaked Ha line but a weak 
H/3 line. The excess emission over the power- 
law around these two lines is likely to be 
contributed by host galaxy light. 

SDSSJIOI54I. 14+594445.2. This is a RL 
quasar exhi biting two pai rs of radio lobes in 
an X-shape ( Cheuneil2007f ) . The origin of the 
X-shape wings in this radio source is unclear. 

SDSSJ101810. 98 +3542 39. 4. This is an 
FSRQ (|Healev et al.ll2007l) . 

SDSSJ102738. 53+605016. 5. The optical 
spectrum of this object exhibits extremely 
broad and double-peaked H/3, and is cata- 
log ed as a dou b le-pea ked emission line quasar 
by IWu fc Liul (|2004[ ). They measured the 
FWHM of H^ of ~ 16, 200 km s"! and a black 
hole mass of log (Mbh/Mq) = 9.649. 

SPSS J 1 03303. 70 +411 6 06. 2. This is an 
FSRQ (|Healev et al.ll2007l) . 

SDSSJ121826. 51+294846. 5. This object, 
also known as Mkn 766 and NGC 4253, is 
a local Seyfert 1.5 galaxy {z = 0.013) and 
is resolved in the UVOT image. Because of 
the strong host galaxy contamination, it is 
difficult to isolate the AGN component. It 
is not included in our UVOT processing list. 
The X-ray spectrum is complicated; it cannot 
be fit by any model we described in Table [21 



23 



This object was previousl y studied in detail in 
the Swift AGN cata l og bvlTurner et alj (|2006l 
2003); iSmmeEaD dioi^^ 

SDSSJ13 5516.54+5612U.7 This is a typ 



ical NLSl (|Zhou et al. 
I999II2OIOI) . 



2006t iGrupe et al 



SDSSJI41927.49+O44513.8. This object 
has a featureless UV spectrum with a few 
narrow absorption lin es. It is a BL Lac ob- 



ject in the catalogs of ICollinge et alJ ( 20051) , 
Plotkin et cdl (|2008l ) and lMassaro et all (|2009l) . 
The narrow absorption lines are most likely 
intervening. 

SDSS J 142921 .87+54 0611.2. This object i s 



class ified as an FSRQ bvlHealey et al.l (|2007| ) , 



and iMassaro et al.l (|2009l ) , and listed as a 
lensed quasar c andid ate by lKing et al.l (|l999i ) . 
Browne et al.l (|2003l ) rejected the lensing hy- 
pothesis based upon surface brightness and 
spectral indices criteria. 

SDSSJ154929.43+023 701.1. This object 
is classified as an FSRQ ('Hcalcy et alll2007l) 
with high polarization (Scarpa &: Falomol 
19971 ). The optical spectrum from SDSS, 
however, contains strong Mg 11 A2798, H/3, 
and even H 7 lines. Even some weak forbid- 
den lines are prominent, e.g., the two [Ne v] 
lines around 3400 A. This is another case in 
which a blazar has a regular broad line quasar 
spectrum. 

SDSSJI6I742.53+322234.3. This is the 
strong ratio quasar 3C 332. The H/3 line is 
relatively weak, but the Ha line exhibits a 
very prominent double peak structure with 
a FWHM of 1 9,600 km s^^ measured by 



Strateva et all (I2OO6I) a nd 23 200 km 



measured by IWu k Lii] (|2004l) . The black 
hole mass is lar ge, log(MBH/-^0) = 9.334 
( Wu fc Liulliool . The Ha line also displays 
long-term profile variability which can be ex- 
plained by a low, smooth, secular c hange in 
disk illumination (jGezari et al.ll2007l ). 

SDSSJl 62901 . 30+400759. 9. This object 



is cla ssified as a bla zar bv IMassaro et ah. 
( 2009h . an FSRQ by iFalcone et al.l (|2004l ). 



and a NLSl galax y bv iBade et al. 



Grupe et al.l (|2004l ): iKomossa et al. 



1995 ) 



20061) . 



The optical spectrum exhibits strong Ha, 
H/3, and a complex of low ionization Fe 11 
emission. 

SDSSJl 7 0231.06+324719.6. T his is a typ- 
ical NLSl (iGrupeet al.l |2004 IZhou et al. 
20061 : iGrupe et al.ll2010l ) 



SDSSJ213638. 58+004154.1. This z = 
1.9414 FSRQ has a strong radio iet dLiu fc Zhang 



2002l:lHealey et al.l2007l:lMassaro et al.l2009[) 
The UV spectrum contains strong Mg 11 A2798, 
C IV A1549, and Cm A1908 emission lines. 
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Fig. 1. — Distribution of our sample in the 
redshift vs. SDSS i band absolute magni- 
tude diagram. The SDSS DR5 quasar cat- 
alog objects are represented by open circles. 
Their distribution is represented by a set of 
linear contours when the density of open cir- 
cles in this two-dimensional space exceeds a 
certain threshold and the plot symbols be- 
gin to overlap. Objects in the Swift quasar 
catalog are represented by red filled squares. 
The lower luminosity limit occurs because 
the SDSS DR5 quasar catalog includes only 
quasars more luminous than Mi = 22.0. 
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Fig. 2. — Four examples of fits to the observed XRT spectra. Intrinsic column densities (A^H.i) are in units of 10^^ cm~^. Spectra 
shown in (b) and (c) are binned as listed in Table [TJ Spectra shown in (a) and (d) are rebinned using the XSPEC command 
setplot rebin (see Section [ 
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Fig. 3.~ Examples of initial SEDs, showing data from Swift UVOT and XRT (black), SDSS (blue) and 2MASS (red). There 
are clear flux offsets between SDSS and UVOT measurements. 
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Fig. 4. — Emission line corrections (EC) of 
the SDSS g band ( Upper panel) and the 
UVOT B band {lower panel) as a function 
of redshift for different filters. Red circles are 
EC performed on real spectra. Blue circles 
are EC performed on the coin posite spectrum 
by IVanden Berk et all (|200ll) shifted to the 
real spectrum redshift. The EC represented 
by blue circles varies with redshift as emission 
lines are shifted within the coverage of a filter. 
These plots indicate that the ECs based on 
real spectra are generally in agreement with 
ECs based on the composite spectrum. 
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Fig. 5. — Distributions of emission line 
correction differences AEC = Edcai — 
ECcomposite in the SDSS g band and the 
UVOT B band based on the results shown 
in Figure m These figures illustrate that we 
can in general obtain consistent ECs based on 
real and composite spectra at different bands. 
The dispersion of their difference is typically 
0.05. 
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Fig. 6. 



Upper: photometric UV slope 



Q!;y,ph vs. spectroscopic UV slope ai^^sp- Low- 
redshift objects (z < 0.8) are in blue and 
objects at higher redshifts (z > 0.8) are in 
black. The dash dotted line is a^/^ph ~ a^^gp- 
The dotted lines represent the median val- 
ues of ckiy,ph = 0.43 and ai/,sp = 0.42, re- 
spectively. Lower: Distribution of differ- 
ences between Q;i,_ph and a,y^sp for low-redshift 
{z < 0.8, in blue) and high-redshift (z > 0.8, 
in black) objects. We also fit these two his- 
tograms with Gaussian profiles. The low red- 
shift sample has a dispersion (a w 0.5) much 
larger than the higher redshift sample disper- 
sion (cr w 0.3). 
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Fig. 7. — Illustration of the photometric shift 
strategy at z = and z — 3. In either case, 
SDSS photometry is interpolated or extrapo- 
lated to the wavelengths of the UVOT filters. 
The SDSS and 2MASS data are then shifted 
by the difference between real and interpo- 
lated luminosities. Because the 2MASS do 
not overlap with the UVOT wavebands, we 
shift the 2MASS photometry with the same 
value as SDSS photometry (see Section [ 
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Fig. 10. — Examples of SED fits with the EXP (solid curves) and TPL (dashed hnes) models. The UV/optical fitting region 
is bounded with a vertical dotted line at 5600 A and a vertical dashed line at Lya. UV/optical data points outside this region 
are not used for SED fitting. Swift UVOT data are plotted in black; SDSS data are plotted in blue; 2MASS data are plotted 
in red. Dotted curves are UV and X-ray components for the EXP model. SDSSJ154929. 43+023701.1 is a case with strong UV 
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Shen et al.l (2008H): red for 



z < 1.9, and blue 
for z > 1.9. We also plot solid, dashed and 
dash-dot lines when the Eddington ratio is 1, 
0.1 and 0.01, respectively. 



40 



80 - 



G 
D 
O 
U 



— I 1 1 r- 



— I 1 r- 




■0.5 



0.0 



EXP model 



I _ _ _ _ _ 



0.5 



1.0 



Fig. 12. — Distributions of Lboi integrated 
from the EXP (thick sohd hne) and TPL (thin 
dashed Hne) models with respect to L^oi cal- 
culated using the BC correction. The black 
shaded region under the EXP model his- 
togram represents objects we flagged as "red" 
which suffer from strong intrinsic absorption 
(see the "Reddened quasars" in Section [O]) . 



l°g Lnjol.model l°g '-'bol.SDSS 



41 



CO 

■l-> 

o 

0) 

IS* 
o 



c 
o 

■l-> 
o 

{-, 

I 
X 



500 


\ ; 


400 


K i 


300 


M ^ 


200 


r \ ; 


100 


""1 

....1 











0.9 


r y-^-'^^^^ ] 


0.8 


■/ ■ 


0.7 


Y \ 


0.6 


r - 


0.5 


' 







20 40 60 80 

Exposure Cut-off (ks) 



100 



Fig. 13. — The clean catalog sample size ( Up- 
per panel) and X-ray detection rate (Lower 
panel) as a function of XRT exposure cut-off 
in kilo-seconds. 



42 



T3 



4 - 



3 - 



2 - 






■ *°; ■ ■ # .°°' ° ° " 



-0.5 



0.0 



0.5 
A(g-i) 



1.0 



1.5 
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(J2003l ). Quasars to the right of the dashed 



line are considered to be dust-reddened. 
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Fiff. 15. — The aox-L ° relation for the 

° 2500 A 

clean sample with a 10 ks XRT exposure cut- 
off ( Upper panel) and the total clean sam- 
ple {Lower panel), including the large clean 
catalog sample (green-|-blue, blue points are 
X-ray detected), and the GIO sample (red). 
Blue points are X-ray detected, while green 
points with arrows are upper limits of unde- 
tected objects. The solid black line is the best 
linear fit to the combined sample using the 
EM method and the solid green line is the 
best linear fit to the cleaned catalog sample 
only. The dot d ashed line is the best fit of 
Just et al.l (J2007I ). The aox in these plots are 
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obtained by fitting photometric data points 
rather than from spectra. The typical error 
bar is displayed at the upper right corner. 
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Fig. 16. — The relation between auv and 
Qx- Objects are selected from the clean cat- 
alog sample (open circles), excluding objects 
with fixed values of a^- The 17 low-redshift 
{z < 0.4) objects are colored in blue and the 
remaining objects with z > 0.4 are in red. 
NLSls are marked with a cross on top of cir- 
cles. We plot green open diamonds on top 
of low luminosity quasars with log L i < 
30.5. For con iparison, we a lso pl ot the AGN 



(|2010[) (open tri- 



sample from iGrupe et al. 
angles). The vertical dotted line marks the 
position where o^x = —1.6. The Spearman 
correlation coefficient for all objects from the 
clean catalog is ps = 0.14 {Pq = 0.125) which 
does not indicate significant correlation. The 
correlation coefficient for low redshift objects 
is Ps = 0.47 (Po = 0.058), for objects with 
Qx > -1.6 it is Ps = 0.22 {Pq = 0.017) and for 
low luminosity quasars it is ps = —0.04 (Pq = 
0.726). 
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Fig. 17. — The relation between ctox and a^ 
for 129 objects (open circles), selected from 
the clean catalog sample, excluding objects 
whose X-ray spectral indices are fixed during 
the fitting process or without L ° mea- 

'^ ^ 2500 A 

sured from UVOT photometry. Data points 
are color-coded depending on redshift with 
z < 0.4 in blue, 0.4 < z < 1.5 in green and 
z > 1.5 in red. Objects with L •>< 30.5 

■' _ 2500 A 

are flagged with larger diamonds in magenta. 
The correlation coefficient of objects in our 
work is ps = -0.06 (Pq = 0.484), which does 
not exhibit a significant correlation (though 
with low confidence level). Th e dotted line 
is the weak correlation found by lYoung et al 



( 2009|) . For comparison, we over plot the GIO 
sample in open triangles. 
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Fig. 18. — The relation between the 

UV/optical spectral index auv and aox for 
217 objects, selected from the clean catalog 
sample, excluding objects whose L2 keV ^^ 
not available or without L = measured 

2500 A 

from UVOT photometry. Data points are 
color-coded in the same way as Fig. [iTl The 
correlation coefficient of objects in our work 
is ps — —0.17, which does not exhibit a sig- 
nificant correlation. Arrows represent quasars 
not detected by XRT. The low redshift (z < 
0.4) quasars in the clean catalog sample are 
color-coded in blue. 
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Fig. 19. — The relation between the Edding- 
ton ratio iboi/iEdd and the UV/optical spec- 
tral index ayv for 247 objects selected from 
the clean catalog sample, excluding objects 
with no bolometric luminosity measurements. 
Data points are color-coded in the same way 
as Fig. [iTl The correlation coefBcient of ob- 
jects in our work is ps = 0.35 (Pq < 10^''). 
The solid straight line is the best linear fit to 
our data by the BCES method. For compari- 
son, we also plot the linear regression results 
of the BLSl, NLSl AGNs and th e comb ine 
of them obtained by IGrupe et al.l ( 20101 ) in 
magenta dash-dotted, blue dash-dotted and 
dotted lines, respectively. 
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Fig. 20. — The relation between Lboi/iEdd 
and ax for 129 objects selected from the clean 
catalog sample, excluding objects with fixed 
values of photon indices. Data points are 
color-coded i n the same w a y as F ig. [TTl The 
AGNs from iGrupe et al.l (|2010l ) are repre- 
sented with black open triangles. The cor- 
relation coefficient of objects in our work is 
Ps = -0.26 (Po < 10"^). The solid black, dot- 
ted and dash-dotted lines represent the linear 
regression r esults to our sainple, th e GIO sam- 
ple and the IShemmer et al.l (J2008I ) sample. 
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Fig. 21. — The relation between the Edding- 
ton ratio iboi/iEdd and ctox for 247 objects 
selected from the clean catalog sample, ex- 
cluding objects whose ctox values are unavail- 
able. Data points are color-coded in the same 



way a s Fig.[T71 The AGNs from lGrupe et al 



( 20101) are represented with black open trian- 
gles. Th e dotted line i s the linear regression 
result bv lGrupeet"aLl (|2010l) . 
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Fig. 22. — X-ray slope F vs. rcdshift (upper panel), L ° (lower left), and L2 ]^q\/ (lower right). In each panel, open eircles 
represent data points in the clean catalog sample. Objects whose X-ray slopes are fixed during the fitting process are excluded. 
Open triangles represent data points in the GIO sample. 



Table 1 
XRT Data Binning Strategy. 



A^xph 


Photon# per Bin 


Statistics 


Flag 


TVxph > 200 


20 


x' 


g(good) 


100 < Wxph < 200 


JVxph/10 


x' 


g(good) 


10 < Wxph < 100 


1 


Cash 


a(acccptable) 


A^xph < 10 


Group Min 1^ 


_2 


u;(wcalc) 


A^xph ~ 






o(Out of FOV) 



In these cases, we group spectral bins with a niininuni of 1 photon per 
bin. 

If the total number of X-ray photons is less than 10, wc do not fit the 
X-ray spectrum but only calculate flux or flux limit. 
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Table 2 
Models used to fit XRT spectra. 



Model 


Qfx 




^H,i 


Objects 


Percentage 


A 


Fixed to 


-1 


Fixed to 


177 


44% 


B 


Free 




Fixed to 


195 


49% 


C 


Free 




Free 


15 


4% 


D 


Fixed to - 


-1 


Free 


12 


3% 
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Table 3 
UVOT Sky Image Flagging Description. 







Source 


Background 


Flag 


Description 


Region 


Region 


-2 


Bad Aspect 


None 


None 


-1 


Out of Image 


None 


None 


1 


Faint 


Default 


Default 





<5 < 0".618 


Default" 


Default" 


2 


0".618 < <5 < 3" 


Customized^ 


Default" 


3 


5 > 3" 


None 


None 


4 


Near edge (ds < 27". 5) 


None 


None 


5 


Not in FOV 


None 


None 


10+ 


27". 5 < dE < 35" 


TBD" 


TBD"^ 


100+ 


dv < 100" 


TBD-^ 


TBD" 



Note. — S is the separation between the SDSS coordinate and 
the Gaussian centroid resulting from fitting the object image (see 
Section I2.3D ; d^ is the distance between the source position and 
the sky image edge; tZv is the distance between the source position 
and the nearest vertex of a sky image. 

^See Section |2.3| for descriptions of the default source and back- 
ground region files. 

Source region circles are centered at the new Gaussian centroid. 

■^Thesc images are passed to the visual inspection process (Sec- 
tion [23} 
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Table 4 
Four types of objects in the raw catalog. 



Type 


UVOT 


XRT 


In Catalog 


BED 


Number 


A 


Y 


Y 


Y 


Y 


637 


B 


Y 


N 


Y 


N 


38 


C 


N 


Y 


Y 


N 


168 


D 


N 


N 


N 


N 


191 


Total 


675 


805 


843 


637 


1034 



55 



Table 5 
Photometric shift strategies. 





UV Bands 


EUV Bands 




Priority^ 


Condition 


Condition 


Mctliod^ 


1 


g+u+U 




E 


2 




g+u+U 


E 


3 


g+B+u 




I 


4 




g+B+u 


I 


5 


r+g+B 




E 


6 


r+V + g 




I 


7 




r+V+g 


I 


8 


V + u+g 




E 


9 


i+r+V 


i+r+V 


E 


10 


g + u+UVWl 




E 


11 




g + u+UVWl 


E 


12 


g+u+UVM2 




E 


13 


g + u+UVW2 




E 


14 




g + u+UVW2 


E 


15 




g + u+UVM2 


E 



The highest priority is 1. 

E for extrapolation; I for interpolation 
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Table 6 
Parameters for black hole mass calculation. 



Rcdshift 



FWHM 



EC 



Reference 



2 < 0.7 


FWHM(H0) 


5100 A 


0.660 


0.53 


0.7 < 2 < 1.9 


FWHM(Mg II) 


3000 A 


0.505 


0.62 


2 > 1.9 


FWHM(C IV) 


1350 A 


0.672 


0.61 



9.26 McLure fc Dunlop (2004*) 

5.15 McLure fc Jarvis (2002) : McLure fc Dunlop (2004) 

3.81 Vestergaard fc Peterson (2006) 



'^Wavelength of continuum monochromatic luminosity. 






Table 7 
Catalog description. 



Column Format 



Symbol 



Description 



1 


A18 


SDSSID 


2 


FIG. 6 


RA 


3 


FIG. 6 


DEC 


4 


F6.4 


z 


5 


F6.2 


Mi 


6 


F6.3 


u 


7 


Al 


Quality Flag 


8 


11 


Catalog Flag 


9 


11 


SED 


10 


11 


A''2MASS 


11 


11 


-^VUVOT 


12 


Al 


XRT Flag 


13 


A6 


QSO Type 


14 


Al 


XRT Model 


15 


12 


Red Flag 


16 


11 


SDSS Fit 


17 


F7.3 


r 


18 


F7.3 


,5-(r) 


19 


F7.3 


<5+(r) 


20 


F7.3 


HR 


21 


F7.3 


log Fob, (G. 3-10 keV) 


22 


F7.3 


logF„„ob=(0.3-lG keV) 


23 


F7.3 


1°S ^1 keV 


24 


F7.3 


<5 log ^2 keV 


25 


F7.3 


logL(0.3-lG keV) 


26 


12 


XRT Detect 


27 


E9.2 


CR 


28 


F8.1 


TxRT 


29 


F5.3 


AfH.G 


30 


F7.3 


AfH.i 


31 


F7.3 


<5-(ArH,i) 


32 


F7.3 


<5+(AfH,i) 


33 


12 


VFlag 


34 


F8.1 


Tv 


35 


F7.3 


log/(V) 


36 


F7.3 


<51og/(V) 


37 


12 


B Flag 


38 


F7.1 


Tb 


39 


F7.3 


log /(B) 


40 


F7.3 


<5 log /(B) 


41 


12 


U Flag 


42 


F7.1 


Tu 


43 


F7.3 


log/(U) 


44 


F7.3 


<51og/(U) 


45 


12 


UVWl Flag 


46 


F8.1 


Fuvwi 


47 


F7.3 


log /(UVWl) 


48 


F7.3 


5 log /(UVWl) 


49 


12 


UVM2 Flag 


50 


F8.1 


TuVM2 


51 


F7.3 


log/(UVM2) 


52 


F7.3 


(51og/(UVM2) 


53 


12 


UVW2 Flag 


54 


F8.1 


ruvw2 



SDSS DR5 Designation hhmmss.ss+ddmmss.s(J2000) 

SDSS right ascension in decimal degrees (J2000) 

SDSS declination in decimal degrees (J2000) 

Redshift from SDSS DR5 quasar catalog 

Absolute magnitude at i band from SDSS DR5 quasar catalog 

BEST PSF u magnitude (not corrected for Galactic extinction) 

Data quality {a/h/c/d') 

1 —in final catalog, —not in final catalog 

1 ^has SED plots, ^no SED plot 

Number of 2MASS photometric points 

Number of UVOT photometric data points 

Quality of XRT data {g/a/w/o) 

Classification of quasar 

Model used to fit XRT spectrum {a/h/c/d) 

1 ^color is "red" , —color is not red 

1 —SDSS spectrum is fit, —SDSS spectrum is not fit 

Photon index between 0.3 and 10 keV 

1(7 lower error bar of photon index 

1(7 upper error bar of photon index 

Hardness ratio 

Observed fiux between 0.3 and 10 keV 

Unabsorbed fiux between 0.3 and 10 keV 

Monochromatic luminosity at 2 keV in erg s~ Hz~ in logarithmic scale 

1(7 uncertainty of monochromatic luminosity at 2 keV 

Integrated luminosity between 0.3 and 10 keV in erg s~ 

1 —detection, —non-detection, —1 —not observed 

Source count rate in 10~ cnt s~ 

Total XRT exposure time in seconds 

Galactic column density in 10 cm~ 

Intrinsic column density in 10 cm~ 

1(7 lower error bar in 10 cm~ of A'^H.i 

1(T upper error bar in 10 cm" of N\i\ 

1 —has V band photometry, —1 —no V band photometry 

Total exposure time (seconds) in V band 

Flux density at Swift V band 

1(T error bar of fiux density at Swift V band 

1 —has B band photometry, —1 —no B band photometry 

Total exposure time (seconds) in B band 

Flux density at Swift B band 

1(T error bar of fiux density at Swift B band 

1 —has U band photometry, —1 —no U band photometry 

Total exposure time (seconds) in U band 

Flux density at Swift U band 

1(7 error bar of fiux density at Swift U band 

1 —has UVWl band photometry, —1 —no UVWl band photometry 

Total exposure time (seconds) in UVWl band 

Flux density at Swift UVWl band^ 

la error bar of fiux density at Swift UVWl band 

1 —has UVM2 band photometry, —1 —no UVM2 band photometry 

Total exposure time (seconds) in UVM2 band 

Flux density at Swift UVM2 band^ 

la error bar of fiux density at Swift UVM2 band 

1 —has UVW2 band photometry, —1 —no UVW2 band photometry 

Total exposure time (seconds) in UVW2 band 
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Table 7-— -Continued 



Column Format 



Symbol 



Description 



55 


F7.3 


log /(UVW2) 


56 


F7.3 


i51og/(UVW2) 


57 


F7.3 


auv.ph 


58 


F7.3 


'5cKUV,ph 


59 


F7.3 


1°SS500A 


BU 


F7.3 


■""g ^2500 A 


bi 


F7.3 


"ox.ph 


62 


F7.3 


<5aox,ph 


63 


F7.3 


logMBH 


64 


F7.3 


logl/bol.EXP 


65 


F7.3 


logibol.TPL 



Flux density at Swift UVW2 band^ 

la error bar of flux density at Swift UVW2 band 

UV spectral index by fitting photometric data 

la error bar of auv.ph 

logL„_j-.^ s by fitting UV photoemtric data in erg s~ Hz 

la error bar of 1 



^2500 A 
calculated using logL„_„„ s and log 



■^2 kcV 



-,ph 



la error bar of a^ 

Black hole mass 

Bolometric luminosity by the EXP model in erg s~ 

Bolometric luminosity by the TPL model in erg s~ 



Note. — Data entry is usually set to —99.9 if unavailable. 

Fluxes are in 10~ erg s~ cm~ A~ , and are corrected for Galactic reddening. 

Calculated using emission line FWHM and corresponding continuum fiux calculated by power-law 
fitting of photometric data points. 

The hardness ratio here is defined as {Nu — Ns) / {Nu + A'^s), in which A''h is the X-ray photon count 
between 1 and 10 keV and A^s is the X-ray photon count between 0.3 and 1 kcV. 
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Table 8 
Data of selected catalog columns for object examples. 



SDSSID 



Quality 
Flag 



XRT XRT 
Flag Model 



log L-2 koV 

[erg s^ Hz^ 



XRT 
Detect 



JVn.i 



flO^ 



flUV,ph 



1°S ^2500 A 
[erg s-i Hz~i] 



000639.20 + 142156.1 


1.3920 


A 


000654.40 + 141442.7 


1.6359 


C 


001141.40 - 004722.6 


1.6483 


A 


001217.08 - 005437.6 


3.6030 


C 


001746.50 - 093546.1 


0.5790 


A 


001904.83 + 003436.5 


2.1178 


A 


001917.31 + 002735.4 


2.4183 


C 


001927.17 + 003539.0 


1.2889 


A 


001927.87 + 003359.9 


1.6234 


A 


001954.60 + 004114.1 


1.9081 


A 


001957.60 + 003936.2 


0.9945 


A 


002303.15 + 011533.6 


0.7285 


A 


002740.38 + 010608.6 


1.5062 


A 


002828.34 -011014.1 


1.1930 


A 


003359.38 + 000230.0 


1.6367 


C 


003409.08 + 000318.4 


1.2015 


A 


003415.77 - 000030.8 


1.9451 


A 


003431.78 - 000957.4 


1.5043 


A 


003435.13 - 000947.8 


1.6711 


A 


003922.44 + 005951.7 


1.9889 


A 


003940.23 + 004241.5 


1.7010 


B 


005446.22 + 140019.0 


0.5015 


A 


005503.52 + 140806.5 


1.6679 


A 


011056.90 + 001912.0 


0.8056 


A 


011119.81 + 002652.0 


1.7475 


A 


011124.42 + 002647.0 


1.0029 


A 


113749.11 + 404913.5 


1.1886 


A 


114502.23 + 595720.0 


1.6385 


A 


145353.56 + 032450.8 


2.4045 


A 


154929.43 + 023701.1 


0.4144 


A 


165004.94 + 313354.6 


1.6948 


A 



9 240"'""'^'^'^ 
z.^^u_Q 220 



2.000 
2.000 



2.000 
2.000 



2.000 



2.000 

1 897+"- 1''^ 

2.000 
2.000 

2.^b8_0.063 

2.000 

1 7iq+0.024 

-L. ( -L»_o,024 

1 911+"'2^i 



27.187 
28.411 
27.050 
28.490 
25.964 
30.320 
27.202 
26.588 
26.936 
26.988 
26.472 
26.807 
26.192 
26.595 
27.870 
26.621 
26.940 
26.823 
26.995 
26.879 



± 0.096 



± 0.062 



±0.152 
±0.114 



± 0.092 
± 0.044 
±0.181 



±0.288 
±0.138 
±0.391 



25.019 ± 
27.164 ± 
26.837 
27.005 ± 
26.402 ± 
26.350 ± 
26.834 ± 
26.965 ± 
27.075 ± 
26.734 ± 



0.089 
0.077 

0.098 
0.105 
0.069 
0.019 
0.128 
0.009 
0.127 



1 


2 K17+1-424 


-0.927 ±0.013 


31.049 ±0.002 


-1.482 ±0.037 





0.000 











0.000 


0.083 ±0.142 


30.086 ± 0.032 


-1.160 





0.000 








1 


0.000 


-0.291 ±0.060 


29.535 ±0.013 


-1.370 ±0.024 





0.000 


-0.290 ±0.073 


30.976 ±0.011 


-0.250 


1 


0.000 








1 


0.000 


-0.673 ±0.089 


30.472 ± 0.013 


-1.491 ±0.044 





0.000 











0.000 


-0.959 ±0.297 


30.652 ± 0.058 


-1.406 


1 


0.000 


-0.513 ±0.074 


29.959 ±0.014 


-1.338 ±0.036 


1 


0.000 


-0.155 ±0.046 


30.430 ±0.010 


-1.390 ±0.017 


1 


0.000 


0.471 ±0.157 


30.021 ±0.031 


-1.469 ±0.070 





0.000 


-0.370 ±0.106 


30.378 ± 0.017 


-1.452 





0.000 








1 


0.000 


-0.039 ±0.051 


30.208 ± 0.008 


-1.377 ±0.111 





0.000 








1 


0.000 


-0.566 ±0.156 


30.444 ± 0.028 


-1.390 ±0.054 





0.000 


-0.338 ±0.155 


30.450 ±0.023 


-1.326 


1 


0.000 


-0.463 ±0.288 


30.723 ± 0.040 


-1.475 ±0.150 


1 




0.037 ±0.040 


30.117 ±0.009 




1 


0.000 


-1.465 ±0.220 


29.577 ±0.040 


-1.750 ±0.038 


1 


647+"-^''2 


-0.872 ± 0.086 


30.626 ± 0.019 


-1.329 ±0.030 





0.000 


-0.412 ±0.016 


30.149 ±0.003 


-1.271 


1 


0.000 


-0.234 ±0.088 


30.724 ±0.013 


-1.427 ±0.038 


1 


0.000 


-0.285 ±0.061 


30.365 ±0.012 


-1.521 ±0.041 


1 


0.000 


0.113 ±0.110 


29.825 ±0.018 


-1.334 ±0.027 


1 


0.000 


-0.111 ±0.045 


30.705 ± 0.009 


-1.486 ±0.008 


1 


0.000 


-0.415 ±0.001 


31.041 ± 0.000 


-1.564 ±0.049 


1 


0.000 


-0.741 ±0.079 


29.948 ±0.015 


-1.103 ±0.007 


1 


^'- ' '^-0.552 


-0.416 ±0.075 


30.268 ±0.016 


-1.356 ±0.049 



Note. — Refer to Tablc[T]for conventions of notations and meanings of flags. Values without error bars are usually because they are fixed (e.g., F and A''H,i) or they represent 
upper limits (e.g., ol^x. and logLo ifov)" Entries without values are because they cannot be measured. The complete catalog content is published in its entirety in the electronic 
version of the Astrophysical Journal Supplement. The portion is shown here for guidance regarding its form and content. 



Table 9 
Quasar types excluded from clean catalog sample. 



Classification Objects 

BAL quasar 50 

BL Lac 1 

Lensed quasar 2 

Extended 1 

RL 97 

Sy 1.8 1 

Dust reddened 50 



Note. — Classifications 

may overlap. 
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Table 10 

Sample selection and properties. 



Sample Catalog GIO Total 



Parent Sample 


637 


88^ 


725 


Small Clean Sample 
X-ray Detection Rate 


214 

85% 


88 
100% 


302 

89% 


Large Clean Sample 
X-ray Detection Rate 


426 
65% 


88 
100% 


514 

71% 






contains 




^The sample in|Grupe 


et al.1 (120101) 


! 92 ob- 



jects, but 4 objects do not have UVOT observations. 
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Table 11 

Correlation and regression analysis. 





Size 


PsiPof 




EM 






BJ 




Sample 


Slope 


Intercept 


Dispersion^ 


Slope 


Intercept 


Dispersion^ 


Small clean catalog 


207 


-0.55 


-0.156 ±0.022 


3.244 ±0.676 


0.122 


-0.157 ±0.022 


3.276 


0.115 


Small clean catalog+GlO 


295 


-0.68 


-0.140 ±0.014 


2.773 ±0.420 


0.113 


-0.142 ±0.013 


2.830 


0.107 


Big clean catalog 


426 


-0.47 


-0.148 ±0.013 


3.018 ±0.398 


0.148 


-0.146 ±0.013 


2.972 


0.140 


Big clean catalog+GlO 


514 


-0.57 


-0.125 ±0.008 


2.338 ±0.246 


0.140 


-0.126 ±0.008 


2.355 


0.131 


J07 


372 


-0.76 


-0.140 ±0.007 


2.704 ±0.212 


0.150 


-0.140 ±0.006 


2.723 


0.132 


J07T* 


289 


-0.66 


-0.134 ±0.010 


2.541 ±0.310 


0.154 


-0.137 ±0.009 


2.612 


0.132 



CO 



Note. — EM: Expectation-Maximization algorithm; BJ: Bucklcy-Jamcs algorithm. 
Standard normal residual. See lLavalley et al.| | |1992|) . 
Kaplan-Meier residual. See lLavalley et al.| | |1992D . 
Spearman rank eorrelation coefficient p^ with confidence level Pq. 
J07 sample after excluding high luminosity quasars with log ^j^ (2500 A) > 31.5. 



